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Preface 


The  use  of  the  Quantitative  Feedback  Technique  (QFT)  in  determining  designs 
that  satisfy  performance  criteria  for  a  set  of  plants  having  structured  parameter 
uncertainties  rests  heavily  on  the  in vertibi lity  of  the  plant  matrix.  While  other  design 
constraints  are  appropriate,  it  is  the  requirement  for  plant  invertability  that  forces 
the  designer  to  either  square-down  or  square-up  a  nonsquare  plant  matrix.  This 
study  investigates  the  use  of  a  frequency  sensitive  weighting  matrix  in  squaring- 
down  a  nonsquare  plant  having  more  inputs  that  outputs. 
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edge,  research  experience,  and  periodic  pep-talks  were  essential  to  the  completion  of 
this  thesis  despite  the  difficulties  that  were  encountered. 

Finally,  I  would  like  to  thank  my  wife  Ingrid  and  daughter  Doris  lor  enduring 
the  sleepless  nights  and  curtailed  levels  of  patience  that  accompanied  this  thesis 
effort. 
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- 

linear  time  invariant. 

mtJ 

- 

element  of  P~l . 

MIMO 

= 

multiple-input  multiple-output. 

MISO 

= 

multiple-input  single-output. 

MP 

= 

minimum-phase. 

NMP 

= 

nonminimum-phase. 

P 

= 

aircraft  roll  rate. 

P'j 

= 

element  ij  of  P(ju>). 

P'; 

= 

element  ij  of  Pe. 

Pv 
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Abstract 

This,  thesis  investigates  the  selection  of  the  frequency  sensitive  weighting  matrix 
needed  to  convert  a  nonsquare  /  x  rn  plant  matrix  into  a  square  plant  matrix  t  hat 
satisfies  Quantitative  Feedback  Technique  (QFT)  multiple-input,  multiple-output 
design  constraints.  The  method  of  pre-specifying  a  desirable  square  plant  matrix  is 
used  in  determining  the  necessary  weighting  matrix.  This  study  assumes  that  the 
nonsquare  plant  matrix  has  more  inputs  than  outputs  and  the  option  to  square-up 
is  not  available  (i.e.,  it  is  assumed  that  inaccessible  states  cannot  be  reconstructed 
due  to  design  limitations). 

The  thesis  examines  several  topics.  The  topics  include  energy  and  power  spec¬ 
tral  densities,  singular  value  decompositions  of  matrices  of  transfer  functions.  QFT 
design  constraints,  QFT  designs  accomplished  on  nonsquare  plants  by  breaking  up 
the  nonsquare  plant  into  a  sum  of  square  plants,  and  both  II2  and  //»- norm  mini¬ 
mization. 

While  it  was  desired  to  pursue  this  study  effort  through  a  comple  e  QFT  design, 
problems  encountered  in  obtaining  the  frequency  sensitive  weighting  matrices  needed 
to  satisfy  a  pre-determined  desirable  plant  matrix  made  this  impossible  in  the  time 
available.  As  such,  neither  the  optimal  squaring-down  compensators  nor  a  final  QFT 
design  is  presented. 

This  study  concludes  with  an  examination  of  problems  encountered  and  rec¬ 
ommendations  for  areas  of  future  study  that  avoid  the  described  problems. 
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WEIGHTING  MATRIX  SELECTION 
FOR  QFT  DESIGNS 


I.  Introduction 

I  his  thesis  focuses  on  the  selection  of  a  frequency  sensitive  weighting  matrix 
Q( j-e)  that  converts  an  original  non-square  /  x  m  plant  transfer  matrix  P(j^)  into 
a  square  plant  matrix  appropriate  for  a  multiple-input,  multiple-output  (MIMO) 
quantitative  feedback  technique  (QFT)  design.  The  plant  is  assumed  to  have  m 
control  inputs  and  l  outputs;  the  number  of  inputs  is  further  assumed  to  exceed  the 
number  of  outputs  (i.e.  states)  available  for  feedback.  The  method  of  pre-specifying 
an  /  x  /  desirable  plant  transfer  matrix  Pe(ju;)  whose  determinant  is  minimum-phase 
(MP)  is  used  in  determining  the  necessary  nonsquare,  m  x  /  weighting  matrix  Q(j~-)- 
r  his  study  assumes  that  inaccessible  states  cannot  be  reconstructed  via  the  use  of 
Luenberger  observers  or  other  estimators  due  to  design  requirements.  The  desired 
Q(j-c)  is  given  as  the  "optimal”  solution  to 


Pe(j^)  =  P(ju?)Q(j*) 


(1.1) 


Synthesis  of  a  QFT  design  for  a  given  nonsquare  plant  transfer  matrix  P(j~')  is 
not.  begun  until  after  the  weighting  matrix  Q(jbu)  is  determined.  Since  the  selection 
process  of  Q(juQ  is  not  an  inherent  part  of  performing  a  QFT  design,  the  use  of  a 
Qfj-G  weighting  matrix  should  increase  design  freedom. 

It  is  anticipated  that  properly  selecting  the  Q(jw)  matrix  yields  a  product 
matrix  P(j*j)Q(ju;)  whose  determinant  is  MP  over  all  plant  parameter  variations 
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an«l  whose  QFT-designed  closed-loop  performance  is  better  than  that  achievable  by 
squaring  down  P(j*cl  without  the  use  of  any  optimization  criteria.  Improvements  are 
expected  because  the  Qf ju.’)  weighting  matrix  acts  as  a  control  input  mixer  to  coor¬ 
dinate  efforts  of  different  control  inputs  in  order  to  increase  the  overall  effectiveness 
of  these  inputs  in  achieving  desirable  plant  responses. 

Different  optimization  criteria  for  use  in  determining  a  Q(j_c)  for  noninvertible 
P(j~)  matrices  ?re  examined.  A  single  QFT  design  problem  is  examined.  While 
the  completion  of  a  QFT  design  was  desired,  the  design  is  not  completed  because  of 
problems  encountered. 

/./  Background 


di(j^)  d2(ju;) 


Figure  1.1.  The  QFT  Compensated  Plano 

figure  1.1  shows  the  QFT  design  matrices  that  are  inserted  as  a  result  of 
synthesizing  a  QFT  design  on  an  /  x  /  plant  matrix.  For  purposes  of  this  study,  all 
system  elements  are  assumed  to  be  expressed  in  the  frequency  domain.  The  system 
elements  are  defined  as  follows: 


u ( j w' )  is  the  control  surface  command  input  as  seen  by  the  aircraft  plant. 
r(jvc’)  is  the  aircraft  command  control  input  provided  by  an  external  source. 
y(jw’)  is  the  aircraft  plant  output. 
di(j.x.’)  and  d2(jw’)  are  disturbance  inputs. 

P(j-0  is  the  uncompensated  plant  matrix. 

F(j  a.’)  is  the  QFT  determined  input  filter  or  prefilter. 

G(jwj)  is  the  QFT  determined  cascade  compensator. 

For  the  QFT  design  problem,  the  output  y(ju;)  is  required  to  track  the  com¬ 
mand  input  r(j^)  and  to  reject  the  disturbances  di(ju;)  and  d2 ( j<-^) .  The  F(ju;)  and 
G(j*y)  elements  are  designed  to  ensure  that  the  sensitivity  of  the  output  to  uncer¬ 
tainties  in  P(j-e)  is  acceptable  and  to  attentuate  the  effects  of  the  two  disturbance 
inputs.  In  addition.  F(ju,’)  is  also  designed  to  achieve  an  acceptable  tracking  re¬ 
sponse.  The  actual  design  and  the  cost  of  feedback  are  closely  related  to  the  extent 
of  plant  uncertainties  and  to  the  narrowness  of  the  performance  tolerances  (7:1.2- 
l.d).  If  a  QFT  design  is  performed  using  di(ju;)  =  0,  then  the  design  problem  is 
referred  to  as  a  type  2  disturbance  problem.  Similarly  if  a  QFT  design  is  performed 
using  d 2 ( )  =  0,  then  the  design  problem  is  referred  to  as  a  type  1  disturbance 
problem. 

Figure  1.2  reflects  the  use  of  a  Q(ju;)  weighting  matrix  with  a  nonsquare  plant 
matrix  P(jw).  A  QFT  design  for  the  MIMO  system  of  Figure  1.2  is  executed  on  the 
product  P(ju;)Q(jw).  As  such,  type  1  disturbance  problems  are  defined  with  di(ju) 
situated  between  Q(j  u>)  and  the  G(ju>)  compensator.  Equation  (1.1)  is  established 
on  the  basis  of  the  open  loop  transfer  function  between  u'(ju;)  and  y(ju>)  which  is 
given  by 


y(jw)  =  PfjwjQCju/Ju'Cjw)  =  Pe(ju>)u'(juj) 
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figure  1.2.  The  QFT  Compensated  Plant  with  Added  Weighting  Matrix 

A  Q(ju-’)  that  minimizes  the  difference  between  the  two  sides  of  Equation  (1.1)  with 
respect  to  some  criteria  is  desired. 

The  selection  of  elements  for  the  Q(ju>)  matrix  is  not  as  straightforward  as 
it  may  seem;  the  selection  process  is  complicated  by  the  existence  of  structured 
uncertainties  in  P(jw).  By  structured  uncertainties,  it  is  meant  that  transfer  function 
elements  of  P(ju;)  vary  so  that  the  plant  is  actually  described  as  a  set  of  plants. 
Conversely,  unstructured  uncertainties  are  modeled  by  only  specifying  boundary 
conditions  on  the  magnitude  of  plant  deviations.  The  set  of  possible  plants  is  denoted 
by  V(jui).  Each  structured  parameter  variation  in  P(jw)  potentially  creates  a  new 
plant  P(ja;)  6  If  the  structured  plant  variations  occur  over  some  continuous 

region,  then  the  resultant  set  of  plants  is  infinite.  The  objective  of  the  QFT  design 
process  is  then  to  satisfy  design  specifications  for  the  set  of  all  plants  within  the 
region  defined  by  the  structured  uncertainties. 

Past  weighting  matrix  designs  have  generally  sought  to  select  elements  of  Q(jtu) 
using  engineering  judgement  and  guesswork  based  on  the  physical  characteristics 
of  the  plant  and  the  importance  of  p,j  elements  of  P(ju)  to  the  achievement  of 
desirable  outputs.  Early  attempts  in  determining  Q(ju>)  concentrated  on  a  Q(ju’)  of 
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constants  only  (1,  16).  Later  attempts  by  Hamilton  and  Phillips  employed  frequency 
dependent  qtJ  elements  in  their  Q(ju:)  (12,  6).  Selection  methods  for  the  transfer 
matrix  elements  of  Q(jw),  however,  remained  a  matter  of  engineering  judgement 
(i.e..  trial  and  error).  For  the  trial  and  error  methods  previously  attempted,  closed- 
loop  performance  results  have  generally  served  to  differentiate  between  effective  and 
ineffective  choices  for  the  matrix  elements  of  Q(ju;). 

Additional  work  is  needed  to  better  define  selection  criteria  for  elements  of 
the  frequency  sensitive  weighting  matrix  Q(jto).  Improvements  in  optimizing  the 
selection  of  Q(>’)  can  also  be  applied  to  other  MIMG  design  techniques  that  use 
weighting  matrices.  An  example  of  such  is  Porter  and  Manganas’  digital  control 
design  technique  that  first  squares-down  a  non-square  plant  matrix  into  an  equivalent 
square  matrix  (13).  A  squaring-down  technique  was  used  by  Phillips  in  his  study 
(12). 

1.2  Problem  Statement 

This  thesis  examines  the  use  of  H2  and  //oo-norm  minimization  as  desirable 
criteria  for  defining  an  optimal  Q(ju;)  compensator  to  Equation  (1.1)  for  a  plant  ma¬ 
trix  P(ju>)  with  structured  uncertainties  subject  to  disturbances  di(jw)  and  d 2 ( ) 
under  the  assumptions  that:  1)  inaccessible  states  cannot  be  reconstructed,  2)  P(ju;) 
has  more  control  inputs  than  outputs  available  for  feedback,  and  3)  that  the  desirable 
plant  matrix  Pe(jw)  is  prespecified.  This  thesis  then  examines  a  design  problem  for 
roll  rate  and  yaw  rate  tracking  controllers  for  a  linearized  model  of  the  AFTI/F16 
aircraft. 

1.3  Review  of  Current  Literature 

Quantitative  feedback  theory  is  essentially  a  frequency  domain  technique  that 
emphasizes  the  use  of  feedback  for  achieving  desired  system  performance  levels  de¬ 
spite  plant  disturbances  and  structured  plant  uncertainties.  Specifics  concerning  the 


1-5 


theoretical  basis  of  QFT  and  problem  applications  are  available  in  (3,  7,  9,  22).  The 
design  technique  is  applicable  to  discrete  and  continuous  systems  for  the  following 
types  of  problems  defined  in  (3:686-695): 

1.  Single-input,  single-output  (SISO)  linear  time-invariant  (LTI)  systems, 

2.  SISO  nonlinear  systems  which  are  converted  to  equivalent  SISO  problems  by 
first  determining  an  LTI  system  that  yields  an  equivalent  set  of  permissible 
outputs  over  the  set  of  all  possible  plants,  disturbances,  and  inputs, 

3.  Multiple-input,  multiple-output  (MIMO)  LTI  systems  which  are  solved  as  a 
set  of  equivalent  multiple-input,  single-output  (MISO)  systems, 

4.  MIMO  nonlinear  systems  which  are  first  converted  to  a  MIMO  LTI  problem 
via  the  same  methods  used  to  convert  SISO  nonlinear  problems  to  SISO  LTI 
problems, 

5.  Distributed  systems. 

The  use  of  QFT  does  not  require  approximations,  linearizations  via  perturba¬ 
tion  models,  or  the  use  of  describing  functions  (9:86). 

In  each  of  the  above  five  types  of  problems,  solutions  to  satisfy  design  param¬ 
eters  are  guaranteed  to  exist  whenever  appropriate  QFT  criteria  are  satisfied.  One 
such  condition  is  that  the  uncertainty  in  P(ju;)  be  known  or  at  least  bounded  (3:695). 
Additional  conditions  apply  and  are  discussed  in  this  study.  It  is  noteworthy  to  in¬ 
dicate  that  QFT  designs  have  been  successfully  applied  to  plants  for  which  a  subset 
of  the  set  of  all  possible  uncompensated  plants  were  unstable  or  nonminimum-phase 
(NMP)  (9,  7). 

For  nonsquare  plants,  the  feedback  control  system  of  Figure  1.2  is  the  simplest 
of  many  possible  schemes.  Horowitz,  et  al,  have  suggested  several  multiple  loop 
structures  that  further  reduce  disturbance  impacts  and  system  bandwidths  in  order 
to  cope  with  large  plant  uncertainties.  However,  while  such  multiple  loop  structures 
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outperform  single-loop  feedback  structures,  the  potential  for  introducing  additional 
noise  sources  and  for  requiring  plant  modifications  to  implement  the  design  are  in¬ 
herent  drawbacks  in  the  use  of  multiple-loop  feedback  structures  (9:31-35). 

Recently,  Yaniv  reexamined  the  MIMO  design  case  (21).  He  suggested  the 
following:  “in  each  design  step,  say  Step  m,  guarantee  that  (1  +  Lm{l))  has  no  [right 
half  plane]  RHP  zeros  for  all  /  6  [  1 , ,  m];  i.e.,  the  same  gm  designed  in  step  m  should 
also  stabilize  the  systems  (1  +  Lm(l)),  (1  +  Lm{2)), ...,  (1  +  Lm(m)).  G  =  diag(glt)  and 
t  =  diag(fit)  will  then  be  the  solution  to  the  problem  stated  ....”  Implementing  the 
change  in  the  design  procedure  ensures  that  all  resultant  diagonal  elements  of  the 
closed-loop  transfer  matrix  have  no  RHP  zeros,  thus  ensuring  that  the  compensated 
system  has  minimum-phase  properties.  The  L ,  loop  transmission  functions  are  the 
L{  =  vuga  of  Chapter  II.  If  the  loop  transmission  functions  have  no  RHP  zeroes, 
then  the  closed-loop  functions  of  the  compensated  plant  are  minimum  phase  (MP). 
Yaniv's  suggestion  applies  to  plants  that  are  basically  non-interacting  (BNIC).  A 
plant  is  BNIC  when  it  is  highly  diagonally  dominant.  In  such  cases,  the  use  of  a 
diagonal  F(ja;)  matrix  in  the  QFT  design  process  is  appropriate. 

1.3.1  Model  Matching  In  order  to  accomplish  the  stated  objectives  -  to  ex¬ 
amine  the  use  of  //2  and  i/oo-norm  minimization  as  desirable  criteria  for  defining 
an  optimal  Q(jcc)  compensator  to  Equation  (1.1)  -  it  is  necessary  to  examine  a 
standard  problem  referred  to  as  the  model-matching  problem.  To  this  end,  the  bi¬ 
lateral  model-matching  problem  is  established  and  then  related  to  Equation  (1.1)  by 
redefining  elements  of  the  model-matching  problem. 
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Figure  1.3.  Model  Matching  Problem 

Figure  1.3  is  the  model-matching  problem  set-up,  where  each  of  the  T:(jw). 
i  €  (1,2,3)  are  fixed  and  Q(ju)  is  an  unknown  compensator  matrix. 

The  objective  of  model-matching  is  to  design  a  stable  controller  Q(ju;)  that 
minimizes  the  difference  between  the  input-output  responses  of  Ti(ju;)  and  the 
cascaded  T2(jw),  Q(jcu),  and  T3(ju;)  transfer  matrices.  Depending  on  the  model¬ 
matching  criterion  used,  the  problem  is  defined  to  either  minimize  the  energy  of  the 
error  z(j u)  for  the  worst  possible  input  v(Jcj)  or  to  minimize  the  largest  error  in  the 
magnitude  of  z(ju>)  for  the  worst  possible  input  v(ju>).  In  the  first  case,  the  problem 
is  termed  an  //2-norm  minimization  problem,  while  in  the  latter  case  it  is  termed 
an  H0 o-norm  minimization  problem.  The  minimization  requirement  is  given  by  the 
following  equation. 


inf  ||$(ju>)||2,oo  =  inf 

QS/Woo  "  ’  QzRHv 


||T1(ju;)-T2(ju;)Q(jW)T3(jw)||2.00 


(1.2) 
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where  inf  is  the  infimum  for  Q(ju)  in  the  space  of  real  rational,  proper  stable 
transfer  functions.  To  be  proper,  the  order  of  the  numerator  polynomial  of  Q(jiu)  is 
less  than  or  equal  to  the  order  of  the  denominator  polynomial.  To  be  real  rational, 
the  transfer  functions  of  Q(ju>)  must  be  expressible  as  a  finite  ratio  of  polynomials 
in  ju>  with  real  constants. 

Solutions  to  the  model-matching  problem  play  an  important  role  in  Hi  and  Hx- 
norm  optimization  (see,  for  example  (5)).  However,  because  this  study  is  directly 
interested  with  the  model-matching  problem,  no  effort  is  made  to  go  into  Hi  and 
H oo-norm  optimization  theory. 

For  this  study  the  model-matching  problem  is  applied  to  the  QFT  problem  of 
determining  an  optimal  Q(jtv)  with  respect  to  Equation  (1.1)  by  redefining  Ti(ju')  = 
Pe(ju>),  T2(jlj)  =  P(jw),  and  Ts(ju;)  =  /.  To  do  so  however,  requires  that  transfer 
function  elements  of  Pe(ju>)  be  stable  (else,  Pe(ju;)  is  in  neither  of  the  two  Hardy 
spaces  of  interest  and  norm-minimization  is  not  applicable).  Also  note  that  for 
structured  uncertainty  problems,  H2- norm  minimization  is  generally  considered  more 
appropriate  than  the  use  of  H^- norm  minimization. 

1.4  Assumptions 

The  following  assumptions  are  made. 

It  is  assumed  that  the  designer  has  either  a  linear  state  space  model  of  the 
system  or  can  establish  an  LTI  equivalent.  A  linear  state  space  model  is  of  the  form 


x  =  Ax  4-  Bu 

(1.3) 

y  —  Cx  +  Du 

(1.4) 
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where  the  A,  B ,  C,  and  D  elements  of  the  equation  are  real  matrices,  and  x  ,  y ,  and 
u  are  vectors. 

An  LTI  equivalent  for  most  nonlinear  plants  can  be  determined  for  the  purposes 
of  QFT  designs.  The  LTI  equivalent  plant  is  basically  established  by  determining 
nonlinear  system  outputs  as  a  function  of  system  inputs  at  particular  frequencies 
and  then  calculating  that  equivalent  set  of  linear  plants  that  produce  comparable 
outputs  under  similar  input  conditions  and  for  the  given  particular  frequencies.  An 
outline  for  determining  LTI  equivalent  plants,  and  conditions  for  their  solvability  is 
given  by  Horowitz,  et  ai,  in  (9:36-44). 

The  AFTI/F16  aircraft  control  example  used  in  this  study  is  developed  under 
the  following  assumptions: 

1.  Aircraft  equations  of  motions  are  linearized  about  the  nominal  equilibrium 
(trim)  condition  of  Mach  0.9  at  20,000  feet.  The  model  is  assumed  valid  for 
small  perturbations  about  the  nominal  conditions.  No  effort  is  made  to 
determine  aircraft  performance  in  regions  outside  the  restrictions 
imposed  by  the  model. 

2.  Mass  of  aircraft  remains  constant. 

3.  Aircraft  is  symmetrical  with  reference  to  a  vertical  plane  aligned  with  the 
longtitudinal  axis  of  the  aircraft. 

4.  Aircraft  thrust  remains  constant. 

5.  Earth  origin  is  taken  as  an  acceptable  inertial  reference. 

6.  Atmosphere  is  a  homogeneous  mixture  at  rest  with  respect  to  the  earth. 

7.  Dimensional  stability  derivatives  are  constant  for  the  given  flight  and  surface 
failure  conditions  modeled. 
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8.  Control  system  failures  can  be  described  by  scaling  the  elements  of  the  B 
matrix.  A  reduction  in  the  magnitude  of  elements  of  the  B  matrix  is  used  to 
represent  gradual  loss  of  control  surface  effectiveness. 

(a)  The  portion  of  the  control  surface  that  has  failed  is  assumed  locked  in  the 
position  associated  with  a  trim  flight  condition. 

(b)  Surface  failures  are  assumed  to  provide  no  other  impact  on  flight  dynamics 
other  than  a  reduction  in  the  effectiveness  of  the  control  surface. 

(c)  Control  surface  failures  are  assummed  to  not  impact  aircraft  stability; 
that  is,  with  control  surface  failures  modeled  by  scaling  columns  of  the  B 
matrix,  the  plant’s  characteristic  equation  does  not  change. 

9.  The  dynamic  behavior  of  the  plant’s  components  in  response  to  inputs  is 
assumed  describable  by  finite-dimensional  Laplace  transfer  functions.  See 
(12,  10). 

10.  Sensor  characteristics  are  such  that,  for  QFT  design  purposes,  the  problem  can 
be  modeled  as  a  type  2  disturbance  problem. 

1.5  Plan  of  Approach 

The  following  is  the  plan  of  approach  for  this  study: 

1.  A  linearized  plant  model  is  selected  for  problem  study.  The  plant  model  has 
three  inputs  (aileron,  rudder,  and  horizontal  tail  deflections)  and  three  out¬ 
puts  of  interest  (sideslip  angle  /?,  roll  rate  p,  and  yaw  rate  r).  The  objective 
of  the  design  is  to  provide  roll  and  yaw  rate  responses  that  satisfy  class  IV 
aircraft  level  1  flying  qualities  (4:412).  For  purposes  of  this  study,  sideslip  is 
not  available  for  feedback. 

2.  Using  the  performance  objectives  specified  in  Chapter  III,  a  desirable  plant 
Pe(ju>)  is  selected  to  satisfy  desirable  plant  characteristics  described  in  (3)  and 


(7)  and  design  criteria  in  (4:412).  The  objective  is  to  determine  a  Q(jw)»  such 
that 

(a)  The  product  of  P(jw)  and  Q(jw)  is  a  square  matrix. 

(b)  The  product  P(jw)Q(jw)  is  non-singular  over  the  entire  range  of  plant 
uncertainties.  Letting  Ve(ju>)  be  the  set  of  all  equivalent  plants  formed 
by  the  product  P(ju>)Q(ju>)  for  any  P(ju>)  G  V(juj),  this  condition  is 
equivalently  stated  by  requiring  that  Ve(juj)  be  nonsingular  over  the  set 
of  possible  equivalent  plants.  Note  that  this  condition  is  applied  to  the 
functions  in  Pe(jcj)  G  Ve{yuj)  and  not  to  evaluations  of  the  functions  at 
particular  frequencies. 

(c)  The  product  P(jw)Q(ju>)  satisfies  a  diagonal  dominance  condition  over 
the  range  of  plant  uncertainties.  This  condition  ensures  that  the  diagonal 
transfer  functions  dominate  the  system’s  response  to  control  inputs  over 
the  entire  range  of  plant  parameter  variations  in  the  high  frequency  range 
(i.e,  above  some  cut-off  frequency).  As  such  the  evaluation  is  conducted 
as  u>  — ►  oo.  A  more  severe  version  of  the  diagonal  dominance  condition 
is  Rosenbrock’s  diagonal  dominance  requirement.  Rosenbrock’s  diagonal 
dominance  requirement  is,  however,  not  a  constraint  that  must  be  satisfied 
prior  to  seeking  a  QFT  design. 

(d)  The  poles  of  P"1^)  lie  in  the  LHP  over  the  range  of  plant  uncertainties. 
In  performing  a  QFT  design,  the  elements  of  P"1^)  are  individually 
inverted.  As  such,  all  poles  of  the  m,j  must  be  in  the  LHP  to  guarantee 
that  the  m~l  are  MP.  To  satisfy  this  criteria  requires  that  factors  of  the 
determinant  de<(Pe(jw))  all  lie  in  the  LHP  over  the  set  of  all  possible 
plants. 

3.  The  Pe(ju>)  chosen  for  use  in  the  optimization  part  of  this  study  is 
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PeCM 


j  Pen  Pe\2 

\  Pej  i  Pei 2 


(1.5) 


with 

peiI  =  562.5/[(jw  +  l0)(—u>2  +  9ju>  +  56.25)] 

Pei2  =  (ju  +  .001)/[(ju;  +  10)2(-u>2  +  9  ju  +  56.25)] 

Pen  =  Pen 

P«2J  =  Pen 

(1.6) 


4.  With  Pe(jw)  specified,  the  Q(ju;)  weighting  matrix  is  determined  in  accordance 
with  optimization  selection  criteria  developed  in  Chapter  II.  The  Q(ju>)  matrix 
is  of  the  format  shown  in  the  following  equation  with  the  q,j  being  the  transfer 
function  elements  of  Q(jiu). 


Q(ju,)  = 


/ 

Ru 
?22 
v  <?32 


qn 
?22 
?32  ) 


(1.7) 


A  QFT  design  is  then  performed  using  the  product  P(ju>)Q(ju>).  Design  objec¬ 
tive  is  to  maintain  class  1  roll  and  yaw  performance  handling  qualities  for  the 
set  of  healthy  and  degraded  plants.  However,  due  to  problems  encountered, 
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this  study  was  not  able  to  obtain  an  acceptable  Q(ja>).  As  such  a  QFT  design 
for  the  plant  described  and  an  evaluation  of  the  design  is  not  accomplished. 


1.6  Scope 

The  method  of  specified  outputs  applies.  The  objective  of  the  Q(ju>)  matrix  is 
to  minimize  the  difference  between  the  original  plant  transfer  matrix  -and  a  prespec¬ 
ified,  more  desirable  equivalent  plant  matrix  Pe(ju;). 

All  analyses  are  performed  in  the  frequency  domain;  digital  considerations  are 
not  addressed.  It  is  noted;  however,  that  both  quantitative  feedback  theory  and  H 2 
and  //oo-norm  minimization  are  applicable  to  both  discrete  and  continuous  functions. 
As  such,  criteria  and  methodologies  examined  in  this  study  are  applicable  to  discrete 
systems. 

1.7  Presentation 

Chapter  II  provides  an  overview  of  pertinent  theoretical  information.  Sections 
of  Chapter  II  address  QFT  theory  while  focusing  on  those  plant  properties  that  must 
be  satisfied  if  a  QFT  design  is  to  be  attempted.  Optimal  means  for  determining 
Q(jw)  are  also  examined. 

Chapter  III  explores  applications  of  QFT  and  weighting  matrix  selection  crite¬ 
ria  to  a  tracking  problem  using  the  lateral-directional  model  of  the  AFTI/F16.  The 
selection  of  Pe(jcj)  and  results  concerning  efforts  to  obtain  Q(ju)  are  presented. 

Chapter  IV  summarizes  conclusions  of  this  study  and  provides  pertinent  rec¬ 


ommendations. 


II.  Theoretical  Background 


2. 1  Introduction 

The  theoretical  background  for  this  study  is  presented  as  follows.  First  energy 
and  power  spectral  density  relationships  both  in  the  time  and  frequency  domains  are 
stated.  Then,  in  order  to  provide  a  basis  for  subsequent  discussions  concerning  sin¬ 
gular  values,  the  singular  value  decomposition  of  matrices  of  constants  and  matrices 
of  functions  in  are  examined.  An  effort  is  made  to  determine  the  availability  of 
transfer  function  descriptions  for  defining  elements  of  the  singular  value  decompo¬ 
sition  of  matrices  of  functions  as  a  means  of  determining  optimal  solutions  to  the 
model-matching  problem. 

With  the  basics  thus  established.  H2  ana  norm  minimization  definitions 
are  made  and  relationships  defined  between  the  norm  minimization  criteria  and 
the  energy  and  power  spectral  densities.  The  Quantitative  Feedback  Technique  is 
next  examined:  both  the  basis  of  the  MISO  equivalent  QFT  design  method  and 
the  source  of  the  QFT  design  criteria  are  established.  Because  of  the  importance 
of  the  minimum-phase  property  to  QFT,  the  impact  of  nonminimum-phase  plants 
and  methods  available  for  determining  when  a  plant  is  NMP  are  next  examined. 
In  addition,  each  of  these  methods  is  examined  in  the  context  of  its  usefulness  for 
establishing  criteria  in  selecting  Q(ju,-)  weighting  matrices  that  guarantee  that  the 
product  P(ju>)Q(j  u.1)  is  MP  for  all  plant  parameter  variations. 

Finally,  each  of  the  QFT  design  criteria  is  reexamined  with  the  goal  of  de¬ 
termining  how  the  criteria  impact  the  selection  of  the  frequency  sensitive  weighting 
mat  rix  Q(ju  ;)  and  how  H2  and  //^,-norms  can  be  used  to  define  a  QCM  'hat  is 
opt  irnal  wit  h  respect  to  Equation  (1.1)  whenever  transfer  function  elements  of  P(ju.’) 
are  known  stable  for  all  plant  parameter  variations. 


2.  2  Spectral  Functions 

Equation  2.1  represents  the  function  of  time  p(t)  as  a  continuous  sum  of  scalar 
exponential  functions  of  frequency  p(ju;);  p(ju;)  is  the  Fourier  transform  or  Fourier 
spectrum  of  p(t).  The  relationship  between  p(t)  and  p(jut)  is  1 


1  r00 

p{t)  =  —  p(jw)e;u"cL,'  (2.1) 

—  7T  J —oo 

/oo 

p(t)e-^dt  (2.2) 

-oo 


A  signal  of  finite  energy  can  be  described  by  a  continuous  spectral-density 
function  by  obtaining  the  Fourier  transform  of  the  signal  (Equation  (2.2)).  The 
Fourier  transform  of  the  signal,  in  turn,  can  be  found  by  setting  s  =  j u.’  in  the 
Laplace  transform  of  p(t),  whenever  both  of  the  following  conditions  are  true. 

1 .  p(t )  =  0  for  t  <  0,  and 

2.  p{t)  has  finite  energy,  that  is,  \p(t)\2dt  <  oo. 

Since  it  is  assumed  that  all  functions  of  time  used  in  this  study  satisfy  the 
above  conditions  ,  the  function  p(jw)  is  obtainable  by  substituting  j  u>  for  s  in  the 
Laplace  transform  of  p{t).  As  a  result,  while  all  functions  in  the  remainder  of  this 
chapter  are  functions  either  of  time  t  or  frequency  the  existence  of  the  Laplace 
transforms  of  p(t)  is  assumed  in  order  to  obtain  p(jw). 

In  cont  rast  to  continuous  functions,  periodic  waveforms  and  the  dc  components 
of  a  signal  have  all  of  their  amplitude  components  at  discrete  frequencies.  To  por¬ 
tray  the  amplitude  components  of  a  periodic  waveform  on  a  spectral  density  graph 

1  The  most  common  convention  in  defining  the  transform  coefficients  is  in  use  here.  See  (U:i  10- 

I  I'-'  ) 
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requires  the  use  of  impulse  functions.  The  area  (weight)  of  an  impulse  is  the  magni¬ 
tude  of  the  time  function's  amplitude  and  thp  spectral  position  of  that  component 
is  the  frequency  at  which  the  impulse  occurs.  (18:83-84). 


If  a  signal  p(t)  represents  a  voltage,  then  p(jca)  has  dimensions  of  voltage 
multiplied  by  time.  The  area  under  the  spectral  density  function  p(ju;),  therefore, 
has  dimensions  of  voltage  because  frequency  has  dimensions  of  inverse  time.  Each 
point  on  a  plot  of  p(j^)  versus  uv  indicates  the  relative  weighting  of  each  frequency 
component.  As  such,  plots  of  p(ju.’)  versus  a;  can  be  considered  as  density  functions. 
The  contribution  of  a  given  frequency  band  to  the  representation  of  p(t)  is  found  bv 
integrating  p(juv)  over  the  bandwidth  area. 

The  square  of  the  absolute  value  of  p(t)  is  called  the  energy  spectral  density  or 
energy  spectrum  of  p(t).  From  the  mathematics  of  the  Fourier  integral,  and  both 
Rayleigh's  and  Parseval's  theorems  (18:85),  the  energy  of  a  signal  E  is  given  by 


e  =  r  i  P(t)\2dt 

J —  Q C 

1  __ 

=  r-  /  P(j-)p(j-)d^ 

2 7T  2-00 

1  f°°  «■» 

=  |p(>’)|2du; 

2  7T  J  —  OO 


(2.3) 


In  this  thesis,  only  real  valued  functions  of  time  are  considered.  Thus,  |p(<)|  = 
p2(t).  By  examining  Equation  (2.3),  it  is  seen  that  the  energy  spectral  density  is 
expressible  in  either  the  time  or  frequency  domain.  Energy  contributions  of  a  given 
frequency  band  are  found  by  integrating  over  the  desired  bandwidth.  The  total 
system  energy  E  is  found  by  integrating  over  (  —  00,00).  Since  the  energy  spectral 
densities  of  real  valued  signals  are  even,  E  is  also  twice  the  value  found  by  integrating 
P’j-’)  over  the  interval  [0,  00). 


In  the  metric  system  of  units,  energy  E  is  measured  in  joules.  Hence,  the 
units  of  the  energy  density  in  the  time  domain  must  be  (joules)2/ second  and  in  the 
frequency  domain  (joules)2/ Hertz. 

Table  2.1  lists  common  energy-storage  elements  that  exist  in  physical  systems, 
and  the  corresponding  energy  equations  for  E.  The  table  provides  insight  into  the 
use  of  spectral  energy  density  functions  in  units  other  than  voltage.  Table  data  is 
extracted  from  (3).  In  the  generalized  sense,  energy  density  expressions  are  given  as 
(quantity)2/  Hertz. 


Table  2.1.  Energy-storage  Elements 


Element 

Energy 

Variable 

Capacitor  C 

Cm 2 

2 

Voltage  v 

Inductor  L 

Li2 

2 

Current  i 

Mass  M 

Mm 2 

2 

Translational  Velocity  v 

Moment  of  Inertia  J 

Ju,2 

2 

Rotational  Velocity  u> 

Spring  K 

h'x2 

2 

Displacement  x 

For  many  applications,  the  power  spectral  density  or  power  spectrum  is  of  more 
interest  than  the  energy  spectrum.  To  describe  differences  between  the  energy  and 
power  spectral  densities  first  requires  an  explanation  of  the  characterictics  of  energy 
and  power  signals. 

An  energy  signal  is  one  that:  (1)  has  finite  energy  content,  (2)  zero  average 
power  (since  the  finite  energy  content  over  infinite  time  yields  zero),  and  whose  form 
is  non-periodic  and  deterministic.  A  power  signal,  on  the  other  hand,  is  one  that: 
( 1 )  has  finite  average  power,  (2)  infinite  energy,  and  (3)  whose  form  is  that  of  either 
a  periodic  or  random  signal  extending  infinitely  over  time. 

Since  energy  and  power  signals  are  inherently  different,  it  follows  that  the 
functions  describing  energy  and  power  verms  time  or  frequency  are  different  for  the 
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two  types  of  signals. 

The  time- averaged  power  of  a  signal  is  given  by  (18:156) 


=  lim  i  r 

T-oo  T 


J=f 


p2(t)dt 


(2.4) 


In  a  manner  analagous  with  the  definition  of  the  energy  E  through  the  use  of 
the  energy  spectral  density,  the  time-averaged  power  is  defined  in  terms  of  a  power 
spectral  density  Sf(ju;)  by 


1  r°° 

2?T  J — oo 


Sf(ju>)dw 


(2.5) 


For  both  random  and  deterministic  signals  whose  power  spectral  densities 
do  not  vary  with  time,  the  power  spectral  density  function  is  obtained  by  taking 
the  Fourier  transform  of  the  autocorrelation  function.  The  autocorrelation  function 
Rj{t)  (18:168-169)  is  given  as 


R}{t) 


1  r  t 

T^of  J=L  P(<  +  T)P(<)^ 


(2.6) 


with  Sf(ju>)  in  Equation  (2.5)  given  as  Sf(ju;)  =  Rf(r)e_ju'Tdr.  For  random, 
stationary  signals  (i.e.,  random  signals  whose  statistics  do  not  change  as  a  function 
of  time),  the  autocorrelation  function  is  given  by 


Rf(r)  =  E{p(t)p(t  +  r)} 


(2.7) 
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where  Rj(t)  is  found  by  time  averaging  (taking  the  expectation  of)  a  statistical 
function  assumed  to  be  ergodic  (  functions  whose  characteristics  are  such  that  time 
averaging  a  single  signal  yields  the  same  statistics  as  obtaining  the  ensemble  average; 
note  that  ergodicity,  in  turn,  implies  stationarity  (but  not  vice  versa))  (18:471-480). 
Certain  nonperiodic  power  signals  are  not  Fourier  transformable;  in  such  cases  the 
power  spectral  density  of  the  signal  is  expressed  in  a  “limiting”  sense  (18:488-4°0). 

Power  spectral  densities  of  periodic  and  random  signals  are  real,  even,  and 
nonnegative  functions  of  frequency-  In  the  metric  system  of  units,  the  power  Pj  is 
measured  in  watts.  Hence,  the  units  of  the  power  density  are  watts/Hertz  for  Sf(jw). 

2. S  Singular  Value  Decompositions  of  Matrices  of  Constants 

1  he  factorization  of  a  real  symmetric  matrix  P  into  P  =  UT,UT ,  UUT  —  /,  and 
with  the  eigenvalues  of  P  in  the  diagonal  matrix  E,  is  known  as  the  spectral  theorem. 
The  matrix  U  is  termed  unitary,  eigenvectors  of  U  are  orthonormal.  The  key  to 
understanding  the  relationship  between  singular  value  decomposition  (which  applies 
to  real  symmetric  matrices)  and  the  spectral  theorem  (which  applies  to  any  matrix 
P)  is  to  recognize  that  the  matrix  products  PPr  and  PT P  are  always  symmetric. 

The  singular  value  decomposition  theorem  provides  that  any  real  /  x  m  matrix 
P  can  be  factored  into  P  =  UT,VT  —  (unitary)(diagonal)(unitary),  where  columns 
of  U  (  /  x  /  )  are  eigenvector.'  of  PPT ,  columns  of  V  (  m  x  m  )  are  eigenvectors 
of  Pr P.  and  where  the  diagonal  entries  of  E  (  /  X  m  )  are  the  real,  non-negative 
values  that  are  termed  the  singular  values  of  P.  For  complex  matrices,  E  is  real  and 
both  0  and  V  are  unitary  (i.e.,  U*U  —  I  and  V'V  =  /  where  (.)*  is  the  conjugate 
transpose  operation)  (17:298,442-443),  Note  that  for  nonsquare  P,  S  is  nonsquare. 

For  example,  consider 
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(  1  +  2;  2  +  3;  3  +  4;  N 

\  4-;'  5  -i  6~i  ; 


Tlie  singular  value  decomposition  for  P  —  UT.V *  is  given  by 


U  = 


'  0.3331  -t-  0.4804; 
y  0.7931  -  0.1533; 


0.1322  +  0.341  iy 
-0.4411  +  0.3911; 


/ 

V*  = 

\ 


0.4190 

-0.8110 

0.4082 


0.5649 

-0.1191 

-0.8165 


0.7108  ^ 

0.5728 

0.4082 


A  connection  between  PP*  and  P*P  is  evident.  The  eigenvalues  of  PP*  and 
P“P  are  the  same.  This  is  readily  seen  by  observing  that  substitution  of  the  sin¬ 
gular  value  decomposition  of  P  into  PP*  and  P“P  results  in  eigenvalue-preserving 
similarity  transformations. 


PP*  =  C/SV*  (PSK*)* 
=  PEE*P* 


(2.8) 


and  similarly, 
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P'P  =  (UEVmy  UEV" 
=  UE*EU* 


(2.9) 


so  that  in  each  case,  similarity  transformations  on  either  EE*  or  E*E  are  obtained. 

Most  applications  involving  singular  value  decompositions  take  advantage  of 
the  numerical  stability  inherent  in  performing  calculations  using  the  decomposed 
form  of  the  P  matrix.  Computational  stability  is  provided  by  the  orthogonality  of 
the  U  and  V  matrices  since  x'U'Ux  =  x'x.  However,  such  a  statement  cannot  be 
made  about  E.  Nevertheless,  E  is  as  accurate  as  possible,  and  it  reveals  the  kind  of 
numerical  stability  problems  which  can  be  encountered  (17:444). 

Another  use  of  the  singular  value  decomposition  of  a  matrix  P  is  in  finding  the 
pseudoinverse  of  P.  The  pseudoinverse  of  a  matrix  provides  the  minimum  length 
projection  solution  to  an  inconsistent  set  of  equations. 

If  the  singular  value  decomposition  of  P  is  given  as  P  =  U EM*,  then  the 
pseudoinverse  of  P,  denoted  as  P+,  is  (17:449) 


P+  =  VE+Um 


(2.10) 


where  the  singular  values  of  P  are  on  the  diagonal  of  E,  and  the  reciprocals,  in  the 
same  order,  are  on  the  diagonal  of  E+. 

Using  the  previous  example  in  this  section,  the  pseudoinverse  of  P  is  given  as 
P+  =  V'E +[/*,  where 
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0.4190 

-0.8110 

0.4082 

0.5649 

-0.1191 

-0.8165 

0.7108 

0.5728 

0.4082 

Vt  _ 


1 


0 


11.0528 
0  ‘ 


0 


.9137 

0 


(  0.3331  -  .4864;  0.7931  +  .1533;  N 

0.7322  -  .3411;  -0.4411  -  .3911;  , 


With  these  results  established,  the  following  claim  is  made.  If  B  is  any  ma¬ 
trix  whose  nonzero  singular  values  are  the  same  as  those  of  P ,  then  there  exists  a 
transformation  matrix  T  that  satisfies  T‘B“BT  =  P*P.  This  claim  parallels  the 
invariance  of  the  characteristic  equation  of  a  square  plant  transfer  matrix  presented 
in  (3:170). 


Let  P  =  UpYjpVp,  and  B  =  Ub'BbVb  such  that  the  nonzero  singular  values  of 
P‘  P  and  B*B  are  equal.  By  assumption,  one  of  the  following  cases  applies: 


case  a.  HpUp 
case  b:  EpEp 


(v*  vD 
0 


case  c:  EpEp  = 


EpEp 


where  the  0  submatrices  are  chosen  to  provide  conformability  between  EpEp  and 
EpE b-  If  case  a,  then 
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SbSb 

VpE*B  SflVp' 

VPZ'BZBVp 

VpVgVBEgEBVgVBVp 

VPVBB’BVBVp 


EpEp 

VpZpUpUpZpVp 

P'P 

P'P 

P’P 


Hence,  select 


T  =  VBVjS 


(2.11) 


The  desired  transformations  for  cases  b  and  c  are  similarly  obtained. 

2.4  Singular  Value  Decompositions  ofP(ju)  Matrices 

Both  MATRIX*  and  MACSYMA  can  be  used  to  obtain  singular  value  de¬ 
compositions  of  matrices  whose  entries  are  constants.  For  example,  the  MATRIX* 
“SVD”  command  is  used  to  obtain  the  singular  value  decomposition  of  any  matrix 
of  constants  (8:4-27). 

MACSYMA,  however,  is  capable  of  performing  numeric  and  symbolic  manipu¬ 
lations.  While  MACSYMA  has  no  command  for  directly  obtaining  the  singular  value 
decomposition  of  a  matrix  of  transfer  functions;  the  “unit  eigenvectors”  command 
is  available  for  obtaining  eigenvectors  and  eigenvalues  of  PPM  and  P* P  needed  to 
construct  the  (/,  E  ,  and  V *  matrices  (19:12-14).  In  using  MACSYMA,  the  E  matrix 
is  obtained  as  the  matrix  formed  by  placing  the  square  root  of  eigenvalues  of  PP* 
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along  the  diagonal.  While  MACSYMA  appears  to  have  the  desired  capability,  it 
must  be  noted  that  for  even  a  simple  2x3  test  case  using  a  matrix  whose  elements 
were  of  the  form  a/(ju>  +  b ),  with  a  and  b  constant,  MACSYMA  requires  several  pro¬ 
cessing  hours  and  returns  a  solution  that  is  extremely  difficult  to  interpret.  For  more 
complex  matrices,  MACSYMA  is  unable  to  return  the  desired  information  because 
of  tne  size  of  the  polynomials  involved  in  performing  necessary  calculations. 

As  such,  if  the  singular  value  decomposition  of  a  matrix  of  transfer  functions 
P(ju>)  is  desired,  it  is  necessary  to  use  numerical  means  to  estimate  elements  of  the 
decomposition. 

The  method  followed  is  to  repeatedly  obtain  the  singular  value  decomposition 
of  a  matrix  P(ju>)  evaluated  at  frequencies  in  an  interval  of  interest.  Points  obtained 
using  the  pointwise  decomposition  of  P(ju>)  are  then  plotted  to  yield  magnitude  and 
phase  plots  of  each  element  in  U(jtu),  V(ju;),  and  X(u;).  The  objective  is  to  then 
estimate  the  transfer  function  elements  of  each  matrix  element  by  using  curve  fitting 
approximations. 

The  following  P(ju>)  is  used  as  a  test  case. 


P(j^) 


(j‘*'+l)(jw+2)(;w+4) 


\ 


0 


_ (jv+ 21 _ 

(ju/+l)(;w+2)(;w+4) 


(jw+2)(jw+3) 

(jw+l)(jui+2)(ju/+4) 


(jw+l)(ju)+2)(ju<+4) 


0U,+4) 

(;w+l)(jw+2)(ju<+4)  / 


The  matrix  P(ju;)  is  evaluated  at  u>  over  the  interval  [  .001,  1000  ]  and  results 
obtained  (see  Figures  2.1  through  2.5). 

The  continuous  lines  of  Figure  2.1  are  plots  of  the  and  a 2(0.')  transfer 

function  elements  of  £(u>)  as  a  function  of  frequency  u>  (the  dotted  lines  are  approx¬ 
imations  to  the  functions).  The  plot  of  <7i(tu)  is  the  upper  plot  of  Figure  2.1  while 
that  of  a 2(0;)  is  the  lower  plot.  Each  of  the  <?\{u)  plots  is  a  nonnegative  function 
with  zero  phase  for  all  u>.  Note  that,  a  real-rational  function  in  u>  cannot  be  used 
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TTTT 


F 


Figure  2.2.  Plot  of  Un(ju;) 
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to  describe  cr1(^')  and  <73(0;).  As  such,  any  description  of  the  a-, {us)  plots  is  accom¬ 
plished  either  by  seperately  describing  the  magnitude  and  phase  of  the  function,  or 
by  using  functions  which  are  not  real-rational. 

The  continuous  lines  of  Figures  2.2  through  2.5  are  plots  of  the  transfer  function 
elements  u;j(j,v).  >  =  1,2,  j  —  1,2,  of  U(jce).  Note  that  uu(ju;)  and  Ui2(j+-’)  cannot 
be  described  as  real  rational  functions  of  ju:.  From  (3:268),  the  phase  difference 
at  =  oc  is  90°'  x  (n  —  re),  where  (n  —  ic)  is  the  difference  in  the  orders  of  the 
denominator  and  numerator  polynomials  of  a  transfer  function.  For  the  given  plots, 
the  phase  plots  of  un  and  of  u  12  indicate  that  (n-  re)  is  odd.  However,  the  zero  slopes 
for  each  of  these  plots  as  — *  oc  indicate  that  the  difference  in  the  orders  of  the 
polynomials  is  zero.  Hence,  any  description  of  the  un(ju;)  and  uio(  ju.’)  must  again 
be  made  either  by  separately  describing  the  magnitude  and  phase  of  the  function, 
or  by  separately  describing  the  real  and  imaginary  parts  of  the  function.  There  can 
be  no  real  rational  function  in  juJ  to  describe  the  first  column  entries  of  U(juz). 

Having  commented  on  the  existence  of  functions  needed  to  describe  £(u.’)  and 
Giju.  ).  a  curve  fitting  approximation  is  obtained  for  each  of  the  provided  plots.  Using 
;UU-)|  to  indicate  the  magnitude  of  the  function,  and  ZU(u,’)  to  indicate  the  phase, 
the  elements  of  S(u;)  and  U(jw)  are  approximated  to  be  as  follows.  Elements  of 
V’fju,’)  (if  desired)  are  similarly  obtained. 

/  r  .  0.994  ll  n  q 

L2 +0.994  J  u  u 

= 

0  f  _0994 — 1 ^  Q 

\  U  L2  +  28.09l  U 

1  I  0  ~02(jw  +  2.~l)(jw+3.75)i  |  0.~ll(ju>+l)(ju/+17)  i 

'  0u/+1.67)(jw+9.8)  '  '  (juv-f-1 .044)(^u>-fl  2.9)  ‘ 

|U(ju;)|  = 

I  0.71 1  ( ;u;  + 1  )(juj  +  1  7)  I  I  0.702Qw  +  2.71  )(jw+3.75)| 

\  1  (ju/+1.044)(jw  +  12.9)  I  1  (ju>+1.67)(jw  +  9.8) 
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The  dotted  lines  of  Figures  2.1  through  2.5  result  from  the  determined  approx¬ 
imations.  In  Figure  2.1,  the  dotted  lines  overlap  the  continuous  (actual)  functions 
quite  well.  If  U*(ju;)  is  required,  it  is  obtainable  as  U*(ju;)  =  |U(ju;)j~l2-  LT(juj). 

Before  closing  this  section,  two  final  notes  are  made.  First,  the  elements  of  the 
E(a,')  and  U(jw)  are  approximations.  It  is  generally  impossible  to  find  a  real-rational 
functions  in  ju>  that  produce  the  curves  obtained  for  elements  of  the  U(jcu),  E(u;), 
and  V(ju;)  matrices. 

Second,  the  functions  used  in  the  test  example  are  simple,  well-behaved  func¬ 
tions.  There  is  no  occurrence  of  a  singular  value  crossing.  If  a  crossing  had  occurred, 
discontinuous  descriptions  of  functions  may  have  been  necessary.  However  since  the 
process  of  curve  fitting  is  imperfect,  function  discontinuities  may  actually  go  unno¬ 
ticed  in  obtaining  a  workable  singular  value  decomposition. 

Hence,  approximations  to  the  singular  value  decomposition  of  a  matrix  of 
transfer  functions  P(jw)  are  obtainable.  However,  because  of  inherent  difficulties 
in  obtaining  acceptable  curve  fitting  approximations  to  polynomials  of  higher  order 
than  are  used  in  the  example  P(jtj),  the  use  of  the  singular  value  decomposition  for 
obtaining  solutions  to  the  model-matching  problem  is  limited. 


2.5  Spectral  Relationships  Using  Singular  Values 

Assume  a  is  a  singular  value  function  of  a  matrix  P,  that  P(ju.’)  is  the  Fourier 
transform  of  P(t)  and  that  P(jw)  is  obtainable  from  P{t).  2 

"‘The  term  P  is  used  in  this  section  to  indicate  that  either  P(f),  or  P(ju>)  is  in  use;  when  there 
is  a  need  to  differentiate,  clarification  is  made. 
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The  function  defining  the  singular  values  <r,  of  P  as  the  positive  square  root  of 
the  eigenvalues  A  of  P*  P  is  given  by 

*,{P)  =  A f  (P'P)  (2.12) 


where  A  satisfies  the  determinant  (det)  equation 

det(\I  -  P'P)  =  0 


(2.13) 


Let  ,4  =  P"P.  Then  for  A  =  [a,j]mxm  and  P  assumed  to  be  /  x  m,  Equation 
(2.13)  is  expanded  to  yield 


( 


(A  —  an) 

+ 

al2 

+  .. 

..  + 

^lm 

a2i 

+ 

(A  —  a22) 

+  .. 

..  + 

«2m 

+ 

®m2 

+  •• 

T  (A  —  a,, 

\ 


-’roXm 


From  the  above  set  of  equations,  the  characteristic  equation  of  A  is  determined.  The 
characteristic  equation  is  given  by 


Am  +  M(m_l)  +  ...  +  A  +  km=0  (2.14) 

Written  in  factored  form,  Equation  (2.14)  becomes 

(A  -  A, )( A  -  A2)  ...  (A  -  Am_i)(A  -  Xm)  =  0  (2.15) 
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The  following  is  known  concerning  the  sums  and  products  of  the  A,  eigenvalues 
of  any  square  matrix  A  (2:208-209): 


m 

y:  A,-  =  trace(A) 

i=i 


(2.16) 


m 

II  ^  =  km 

1=1 


(2.17) 


The  product  of  two  matrices  Mi  and  M2,  Mi  =  [ar,J]n*m  and  A/2  = 


is  defined  by 


K-]nxm  •  l(3,}]m*q  =  hv]nxq 


(2.18) 


where 


E 


fc=i 


atikPkj  =  lij ;  i  =  1,2,...  n,  j  =  1,2,...  q 


(2.19) 


For  A  =  P* P  ,  P  =  \pij]qxn  and  P *  =  [p;i]nX?,  Equation  (2.19)  yields 


<7 

7.J  =  1,7  =  1,2, ...  n 

fc= l 


(2.20) 
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To  obtain  elements  of  the  trace(A)  =  trace(Pm P),  set  i  =  j  and  obtain 


<? 

7..  =  E  PkiPkx  *  =  l,2,...,n 

/t=i 


(2.21) 


so  that  trace(A)  =  trace(P* P)  is  given  by 


n  n  g 

trace(A)  =  E  7 a  =  E  E  PkiPki 


i  =  1  i=l  fc=l 


(2.22) 


Equations  (2.16)  and  (2.22)  are  combined  to  yield 

E  A<  =  E  E  PkiPki  (2.23) 

i=l  i=l  *= 1 

Hence,  if  the  ph ,  are  replaced  by  their  Fourier  transforms,  it  can  be  concluded 
that  the  sum  of  the  squared  singular  values  of  P(ju>)  (the  matrix  whose  elements  are 
the  Pki(jw)  )  evaluated  at  any  frequency  equals  the  system’s  total  spectral  energy 
density  function  (the  sum  of  the  PkiCjkOPkiO^)  evaluated  at  the  same  frequency). 

Integrating  Equation  (2.23),  and  comparing  it  to  Equation  (2.3)  indicates  that 
the  total  energy  E  in  a  system  is  given  by 


E 


1  too  n_  ^ 

rj  EZAm 

ut  J-O o  )=1 fc=1 


1  roo  n 


(2.24) 
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1  r°° 

—  /  trace(A(t))dt 
2  7T  7  —  00 


(2.25) 


For  fun-tions  in  the  frequency  domain,  E  is  given  by  3 


E 


1  roo  "  9 

FT/  EEPki(j^)Pki(jw)d- 
27T  y_oo  ,  =  1  fc=1 

1  r°°  n 


(2.26) 


Spectral  Relationships  to  Norm  Minimization 

To  begin  this  section,  the  Hardy  spaces  known  as  Hi  and  Hoo  are  defined. 

The  H?  space  consists  of  all  compex-valued  functions  P(ju>)  that  are  analytic 
in  the  open  right  half  plane  and  that  satisfy  the  following  condition 

(2-  J  frace[?*(jiu)P*(ju;)]du^  <  oo  (2.27) 


The  left  hand  side  of  the  inequality  is  defined  to  be  the  F^-norm  of  P(ja>),  and 
is  denoted  as  || P(jov) jj 2- 

Alternately,  H0 0  consists  of  all  complex-valued  functions  P(jw)  which  are  an¬ 
alytic  and  bounded  in  the  open  right  half  plane.  By  bounded,  it  is  meant  that  there 
exists  a  real  number  6  such  that 

Eigenvalues  of  P‘  P  equal  the  singular  values  of  P'P  (proof  given  later).  Hence,  the  Aj(jw)  in 
Equation  (2.26)  are  nonnegative  functions  of  jui. 
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Vmax [P(jw)]  <  b 


(2.28) 


where  amax  represents  the  largest  singular  value  of  P(ju>).  The  least  such  bound  b  is 
the  H0 o-norm  of  P(ju;),  denoted  ||P(ju;)||oc>.  Equivalently 

l|P(jw)IL  =  supcrmax[P(ju;)]  (2.29) 


where  sup  over  l j  is  the  supremum  of  the  maximum  singular  value  crmax  of  P(jw) 
over  all  frequencies. 

Relationships  between  the  H2  and  H^- norms  and  the  power  and  spectral  den¬ 
sities  are  now  defined.  Equations  (2.26)  and  (2.27)  are  combined  to  obtain 


1  f°°  71  _  9  ^ 

l|P(j")lla  =  (tt  f  £  £pid(j^)Pki(j^)<M2' 

2tt  7-00  i=1  k=x 
1  n  1 

= 

=  Ve 


(2.30) 


Hence,  for  any  signal  v(jw),  the  //2-norm  of  v(jw)  is  interpretable  as  the  square 
root  of  the  energy  of  the  signal.  This  fact  is  noted  by  Francis  (5:6).  That  the  //2- 
norm  is  also  related  to  the  trace  of  the  power  spectrum  is  obvious  after  substituting 
Equation  (2.12)  into  Equation  (2.29). 

Now  the  standard  model-matching  problem  is  considered.  The  objective  of 
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Figure  2.6.  Model-Matching  Problem 

model-matching  is  to  minimize  ||$(jw)j|  where  ||$(ju>)||  =  ||Ti(ju;)  — T2(ju;)Q(ja;T3(jw)||, 
Ti(a.’)  fixed  and  Q(ju;)  unknown  (see  Figure  2.6)  .  By  comparing  the  model-matching 
problem  with  Equation  (2.30),  it  is  concluded  that  //2- norm  minimization  minimizes 
the  total  energy  of  $(ju).  A  similar  relationship  exists  for  the  //oo-norm. 

Squaring  both  sides  of  Equation  (2.29)  and  letting  cr^ax  denote  the  square  of 
the  maximum  singular,  equations  (2.12)  and  (2.29)  are  applied  to  obtain  the  following 
series  of  equalities. 


IIPCMIII 


(sup  CTmax  [P(ju>)])2 

Uf 

SUP  °max  IP(jw)] 

U / 

sup  Xmax  [P*(ju;)P(jw)] 

U / 

sup  xL  [P*(jw)P(jw)]*  [P-(jo/)P(jw)j 

Ul 

sup  (7max  [P*(jw)P(ju;)] 

Ul 


(2.31) 
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In  addition,  based  upon  equalities  between  intermediate  steps  in  Equation  (2.31),  it 
is  also  concluded  that 


sup  amax  [P*(ju)P(jLj)]  =  sup  AmaI  [P*(ju;)P(ju>)]  (2.32) 


Equation  (2.32)  indicates  that  the  singular  values  of  any  Hermetian  matrix  A  = 
[P'(jcj)P(ju)]  equal  the  eigenvalues  of  A.  As  such,  the  eigenvalues  of  [P*(ju)P(jt<;)] 
are  nonnegative. 

Next,  the  product  of  a  unitary  matrix  U(ju;)  and  a  matrix  P(jcv)  is  formed. 
The  //0 Q-norm  of  the  product  is  given  by 


l|u(jw)  pomil 


sup  [U(ju)  P(jul)] 

U> 

sup  aLx  [U(ju;)  P(ju,)]*  [U(jw)  P(jw)] 

u/ 

sup  \LX  [P*(ju>)  U*(jw)U(jw)  P(jw)] 

C Jj 

sup  \Zxax  [P*(j^)P(jw)] 

Uf 

i|P(j^)!!oo 


Similarly,  it  can  be  shown  that  ||P(ju>)U(ju>)||00  =  ||P(jw)||00. 

As  a  result,  it  is  concluded  that  the  H0 0-norm  of  matrices  formed  by  the  pre- 
or  post-multiplication  of  any  matrix  P(ju>)  by  unitary  matrices  is  the  f/oo-norm  of 
P(jw)  alone. 

A  parallel  result  concerning  the  //2-norm  of  unitary  matrices  can  also  be  shown. 
That  is, 
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l|U(ju>)  P(ju)||, 


(i-  r  tract\{V( jwJPfju)]- [U(ju,)P(jo,)]]^)i 

Z7T  7  —  co 

1  7°°  l 

(—  /  (raa[p-(ju>)U-(jw)  C(ju»)P(ju,)]<iw)> 

Z7T  7  —  oo 

1  /*oo  j 

(—  /  trace[P*(jw)  P(jw)]dw)a 

Z7 T  7  —  00 

l|P(j«)ll2 


and, 


||P(jw)  U(jw)||2 


1  r°°  i 

(-  /  (race[(P(jw)U(ju.)r[P(ju.)U(jw)])<^)‘ 

Z  7T  7—oo 

(-L  /”  /race[[P(jw)U(ju))]  [P(jw)U(juOn«M’ 

Z7T  7  —  oo 

l|P(jw)||a 


where  the  fact  that  trace(AB)  =  trace(BA)  is  used.  Thus  unitary  matrices  are  both 
//2  and  H0 0-norm  invariant.  Multiplying  any  matrix  by  a  unitary  matrix  does  not 
change  either  its  i/2  or  //co-norm.  In  addition,  it  is  also  concluded  that  inherent 
differences  between  energy  and  power  signals  makes  the  choice  of  whether  to  use  //2 
or  Hw norm  minimization  criteria  application  dependent. 

2. 7  Boundaries  on  Model  Matching  Problem  Solutions 


A  sufficient  condition  for  the  existence  of  a  solution  to  the  model  matching 
problem  is  provided  by  Theorem  1.1  (5).  The  theorem  states  that  an  optimal  Q(ju>) 
to  the  model  matching  problem  exists  if  the  ranks  of  the  T2(jw)  and  TaQw)  are 


constant  for  all  0  <  ui  <  oo. 


Now  again  consider  Equation  (2.1),  the  minimization  equation  using  the  //2  or 
Hoc- norm.  As  both  the  //2  and  the  Hoc-norm  are  normed  spaces,  they  satisfy  the 
Schwartz-inequality  and  hence 

11^(^)112,00  =  mow)  -  T2(jW)Q(ju/)T3(jw)lla.oo  (2.33) 

implies  that 

II ^ ( J^) II 2,oo  <  ||Ti(ju;)||2i00 ||T2(ju;)Q(ju’)T3(ju;)||2i0O  (2.34) 

Using  the  above  equation  and  inequality,  we  are  able  to  define  both  a  lower  and  an 
upper  bound  for  model  matching  problem  solutions.  The  lower  bound  is  given  when 
Ti(jw)  =  [T2(jw)Q(ju;)T3(ju;)]  in  which  case  ||$(ju>)||2iCO  =  0. 

The  upper  bound  is  given  by  setting  Q(jcu)  =  0  in  which  case  || ^ ( jcti) ||oo  = 
||Ti(ju;)||2i00.  Hence,  if  a  solution  exists,  then  the  optimal  solution  can  be  sought  via 
an  iterative  process  on  selections  of  Q(ju>)  £  RHoo ,  since  the  following  inequality  is 
always  satisfied 

0  <  ${juOlkoo  <  ||T1(jW)||2i00  (2.35) 

The  Quantitative  Feedback  Technique  is  now  examined  before  relating  QFT 
to  the  use  of  //2  and  H ^  minimization. 
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3.S  The  Quantitative  Feedback  Technique 

Figure  2.7  represents  an  /  x  /  MIMO  closed-loop  system  in  which  F(ju>)  is  the 
matrix  of  prefilter  transfer  functions,  G(jo>)  is  the  matrix  of  compensating  transfer 
functions,  and  Pe(jo;)  =  P(jw)Q(jw)-  Note  that  the  feedback  matrix  is  set  to  /,  since 
an  equivalent  cascade  compensator  with  prefilter  is  determinable  for  any  nonunity 
expression  that  could  be  used  here. 


Figure  2.7.  The  QFT  Compensated  Equivalent  Plant 

When  performing  a  QFT  design,  the  compensator  G(jcj)  is  typically  assumed 
diagonal.  The  use  of  nondiagonal  G(jw)  matrices  has  been  used  in  the  past  to  give 
designers  more  flexibility  ir^ their  designs  (3:698).  However,  the  use  of  nondiagonal 
G(jwv)  increases  the  work  associated  with  determining  G(jcj). 
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From  Figure  2.7,  the  following  equations  can  be  written: 


y(jw) 

=  Pe(j^)C(K)e(ju;) 

(2.36) 

e(>’) 

=  F(jw)r(ju;)  -  y(ju>) 

(2.37) 

y(j-) 

=  [/  +  Pe(>-)G(ju,)]-1Pe(j^)G(»F(ju-)r(j.v) 

(2.38) 

Equation  (2.38)  above  gives  the  MIMO  system  control  ratio  T(ju,')  relating 
r(jwuf)  to  y(ju>)  as 

T(jw)  =  [/  +  Pe(ju;)G(j^)]-1Pe(ju;)G(jW)F(ju;)  (2.39) 


The  QFT  design  procedure  requires  finding  the  unknown  transfer  function  el¬ 
ements  of  G(jtu)  and  F(jw),  by  using  performance  specifications  to  define  tolerance 
bounds  on  the  transfer  function  elements  of  [I  +  Pe(ju;)G(ju>)]_1Pe(ju;)G(ju;)F(ju;). 
To  do  so,  a  mapping  between  a  MIMO  plant  and  a  set  of  equivalent  multiple-input, 
single-output  (MISO)  control  systems  is  defined.  The  mapping  results  in  / 2  equiva¬ 
lent  systems  for  a  MIMO  plant  Pe(jiu)  assumed  to  be  l  x  /,  with  each  transfer  function 
of  the  equivalent  MISO  system  assumed  to  have  two  inputs  and  one  output.  Of  the 
two  inputs,  one  is  designated  the  desirable  input,  while  the  other  is  designated  the 
“disturbance-1  input  (3:696-699). 

The  justification  of  the  MISO  equivalent  method  follows  (7)  and  (3).  The 
in\erse  of  a  square,  (assumed  invertible)  plant  matrix  is  represented  by 
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mn  m,  2 

m  21  m2  2 

^ml  ^mm  f 

where  the  mu  are  elements  of  the  adjoint  of  Pe,  divided  by  det(Pe).  The  /  x  /  M1S0 
equivalent  plant  is  the  matrix  formed  by  inverting  individual  elements  of  P ~ 1  f j ) 
so  that  each  element  of  V7  satisfies  vt]  =  (m,j)_l.  The  MISO  equivalent  matrix  is 
given  by 


mlm 


'  (mu)  1  (m12)  1 

(m2i)_l  (m22)“l  • 

•  •  (mim)'1  ^ 

_ 

'•’ll  t’12 

v2\  v22 

'  •  ^lm  ^ 

^  (^mi)'1  (mm 2)_1  . 

( rnmm )  J 

\  "ml  Vm2 

‘ 

^mm  J 

The  matrix  P“1(jin)  is  next  partitioned  into  two  parts;  the  first  part  D(ju.’) 
contains  the  diagonal  elements  of  and  the  second  part  B(jw)  contains  the 

balance  of  Pg(ju;).  That  is 


=  D(ju;)  +  B(ju;) 

=  [mul 


(2.40) 


Both  sides  of  Equation  (2.40)  are  next  premultiplied  by  [I  +  Pe(ju>)G(ju>)],  to  obtain 


[/  +  Pe(ju-)G(ju;)]T(ju,)  =  Pe(jn.’)G(jtn)F(jcu)  (2.41) 
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Next,  both  sides  of  Equation  (2.41)  are  premultiplied  by  P~1(ju.’)  (assumed  invert¬ 
ible),  to  obtain 


(PeV)  +  G(ja>)]T(ju;)  =  G(jw)F(ju;)  (2.42) 


1  he  equality  Pe^ju,’)  =  D(ju ;)  -f  B(ju>)  is  now  substituted  into  Equation  (2.42)  to 
yield 


[[D(jw)  +  B(jw)]  +  G(jo;)]T(ju>)  =  G(ju,)F(jw)  (2.43) 


Multiply  out  the  terms  on  the  left  side  of  the  equation.  Subtract  B(ju?)T(jw) 
from  both  sides  and  then  premultiply  by  (D(jw)  +  G(ju;)j_I,  which  is  invertible  for 
G(ju>)  assumed  to  be  diagonal.  Thus 


T( jw)  =  [D(jw)  +  G(ju.’)]_1[G(ju;)F(ju;)  -  B(jw)T(ju>)]  (2.44) 


This  last  equation  is  used  to  define  the  desired  mapping.  The  proof  of  the  fact 
that  designs  based  on  MISO  equivalent  systems  yield  satisfactory  MIMO  designs 
for  the  original  plant  is  based  on  Schauder’s  fixed  point  theorem.  The  theorem  is 
invoked  by  defining  a  mapping  Y(T)  -  [D  +  G\~'[GF  -  BT\  where  each  member 
of  T  is  from  an  acceptable  set  of  control  ratios  satisfying  the  design  requirements 
(3:099-700). 
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Individual  elements  of  [D  +  G]  1  are  scalars  multiplying  rows  of  [GF  —  DT ], 


with 


[d,i  +  ga\  1 


[-L+9,.]-' 

vg 

vg 

1  +  vg  gu 


(2.45) 


For  a  unit  impulse  forcing  function,  the  yt]  elements  of  the  mapping  Y(T)  take 
on  the  form 


7— 7— -[S  iifij-'EhijUj] 

I  ~r  diiQii 


Vii 


,  [ggf\j  vijUj] 

1  +  vugu 


(2.46) 


Equation  (2.46)  is  also  the  control  ratio  of  the  ith  MISO  system.  The  transfer 
function  given  by 


vg 

1  +  V„g, 


(2.47) 


represents  the  "desired”  ith  output  with  respect  to  the  desired  jth  input  while  the 
transfer  function 


2-32 


1  +  vaga 


(2.48) 


relates  the  ith  output  to  the  jth  '‘disturbance”  input.  The  objective  of  the  design  is 
then  to  have  each  loop  track  the  desired  input  while  minimizing  the  outputs  due  to 
the  disturbance  inputs. 

For  plants  that  are  basically  noninteracting  (BNIC)  (i.e.,  highly  diagonally 
dominant),  the  magnitude  of  functions  represented  by  Equation  (2.48)  is  smaU  rel¬ 
ative  to  those  represented  by  Equation  (2.47).  For  such  systems,  QFT  designs  can 
be  pursued  using  diagonal  F(ju;)  prefilters. 

For  disturbance-rejection  problems,  the  ytJ  responses  are  designed  to  be  less 
than  some  bound  k\  the  boundary  value  with  respect  to  the  ith  input  and  the  jth 
output,  kij,  is  then  used  to  obtain  the  desired  MIMO  design  by  obtaining  transfer 
functions  that  satify  k\:  >  1 J [ZTfc#.  vX}ttl\.  Solutions  are  obtained  by  reorga¬ 
nizing  the  inequality  into  a  form  that  yields  the  loop  transfer  function,  and  then 
following  the  QFT  process  given  in  (7). 

It  has  been  shown  that  a  realistic  definition  of  optimum  for  purposes  of  QFT 
designs  on  LTI  systems  is  the  minimization  of  the  high-frequency  loop  gain.  It 
has  also  been  proven  that  this  optimum  lies  on  the  composite  boundary  formed  by 
determining  the  maximum  of  two  error  sources:  plant  deviations  due  to  disturbance 
inputs  and  tracking  errors.  It  has  also  been  proven  that  the  optimum  is  unique 
(7:2-26).  The  optimal  design  is,  however,  not  achievable.  Instead,  designs  that  yield 
responses  that  are  relatively  close  to  the  composite  boundary  without  crossing  it  are 
deemed  acceptable,  although  they  encompass  some  amount  of  overdesign. 

The  QFT  process  is  simplified  if  the  mtJ  of  P~l  are  MP.  This  is,  however,  not 
a  QFT  requirement,  since  the  vtJ  =  (m.y)-1  need  not  be  stable  for  development  of  a 
QFT  design. 
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However,  since  poles  of  the  mtJ  of  P~l  appear  as  zeros  of  the  vtJ  functions  used 
in  the  QFT  design,  the  use  of  QFT  in  its  simplest  form  requires  that  all  poles  of  the 
niij  be  in  the  LHP.  Hence,  it  is  important  to  consider  the  impact  of  not  having  all 
LHP  zeros  on  a  QFT  design  before  proceeding  further. 

2.8.1  The  Impact  of  NMP  vtJ  on  QFT  Designs  A  vtJ  transfer  function  of  a 
MISO  equivalent  plant  is  NMP  if  a  zero  of  the  transfer  function  lies  in  the  RHP. 
Systems  with  RHP  zeros  cause  considerable  difficulty  in  control  designs;  the  zeros 
have  a  marked  influence  on  the  nature  of  a  system’s  time  response.  Whenever  there 
is  a  RHP  zero,  the  initial  time  response  of  the  plant  is  negative,  even  though  the 
steady  state  value  is  positive.  In  addition,  branches  of  the  root  loci  that  start  at 
open-loop  poles  must  end  at  the  zeros  or  oo.  For  any  simple  feedback  system  that 
may  be  implemented,  the  closed-loop  poles  lie  somewhere  on  the  root  locus  between 
a  starting  pole  and  either  an  ending  zero  or  oo.  The  location  of  the  closed-loop 
poles  depends  on  the  selected  value  of  feedback  gain.  As  the  selected  value  of  gain 
is  increased,  poles  of  the  closed-loop  system  approach  the  ju>-  axis.  If  the  gain  is 
selected  large  enough,  closed-loop  poles  in  the  RHP  are  obtained,  and  the  closed-loop 
system  becomes  unstable. 

For  NMP  V{j,  tolerances  on  both  the  phase  and  the  magnitude  of  the  tracking 
response  must  be  specified;  for  MP  u,;,  only  the  magnitude  of  the  tracking  response 
is  specified.  While  it  is  considerably  easier  to  work  with  MP  systems,  it  is  not 
an  absolute  requirement  in  order  to  perform  a  QFT  design  (3:692).  QFT  designs 
are  typically  pursued  only  for  MP  systems,  since  it  is  for  such  systems  that  QFT 
guarantees  a  solution  meeting  stated  performance  objectives  for  the  entire  set  of 
possible  plants  (assuming  tolerances  have  not  been  too  tightly  specified). 

The  occurence  of  NMP  behavior  is  a  function  of  where  actuators  and  sensors 
are  mounted.  System  zero  locations  are  changed  if  the  location  of  the  actuator  and 
sensor  devices  is  changed.  Hence,  the  NMP  nature  of  a  vt]  can  be  addressed  as  a  need 
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to  alter  the  systems’s  configuration  by  relocating  control  devices.  However,  once  the 
location  of  a  set  of  actuator  and  sensor  devices  is  decided  upon,  any  NMP  behavior 
and  inherent  control  design  difficulties  that  result  must  be  dealt  with  accordingly 
(15:2-36). 

2.8.2  Methods  for  Determining  if  a  Matrix  is  NMP  The  current  literature  of¬ 
fers  a  number  of  ways  to  determine  if  any  are  MP.  Three  methods  are  examined. 

2.8.2. 1  Routh's  Stability  Criterion  Routh’s  criterion  is  a  simple  method 
for  determining  the  presence  of  RHP  roots  without  actually  solving  for  the  roots. 
The  criterion  states  that  the  number  of  roots  of  the  polynominal  with  positive  real 
parts  is  equal  to  the  number  of  sign  changes  of  array  coefficients  in  the  first  column 
(3:185-191). 

While  Routh’s  criterion  is  primarily  used  for  determining  system  stability,  it 
can  be  used  to  determine  if  any  have  zeros  in  the  RHP  by  applying  the  criterion 
to  the  denominator  polynomial  of  def(Pe(ju;)).  The  application  of  Routh’s  criterion 
is  straightforward,  when  Q(ju>)  is  known  and  P(jw)  is  fixed.  However,  if  the  plant 
varies,  and  if  Q(jw)  is  unknown,  then  the  application  of  Routh’s  criterion  becomes 
increasingly  difficult. 

2. 8. 2. 2  Evaluation  Based  on  Step  Response  The  author  of  (20)  offers 
a  simple  necessary  and  sufficient  condition  for  a  stable  SISO  plant  to  exhibit  an 
undershooting  step  response.  Specifically,  the  author  shows  that  undershoot  occurs 
if  and  only  if  the  transfer  function  of  a  stable  SISO  plant  has  an  odd  number  of  real 
RHP  zeros. 

Consider  a  system  with  a  strictly  proper,  real-rational  transfer  function  p(ju/). 
Further,  suppose  the  system  is  stable  and  has  a  steady  state  value  not  equal  to  zero, 
bet  y(t)  denote  the  time  response  of  the  system  in  response  to  a  step  input.  Then 
for  an  assumed  stable  system,  y(oo)  is  well-defined  and  equals  p(ju;)  evaluated  at 


u>  =  0. 

The  author  show  that  undershoot  in  stable  systems  in  response  to  a  step  input 
occurs  only  if  the  steady-state  value  has  a  sign  opposite  from  that  of  its  first  nonzero 
derivative.  To  define  when  undershoot  occurs  first  requires  defining  the  relative 
order  of  a  polynomial.  The  relative  order  of  a  transfer  function,  denoted  as  r,  is  the 
difference  between  the  orders  of  the  numerator  and  denominator  polynomials  of  the 
function.  In  (20),  a  plant  is  shown  to  exhibit  undershoot  only  if  the  product  of  the 
rtk  derivative  and  the  steady  state  value,  yr(0)  3/(00)  is  less  than  zero. 

The  evaluation  method  described  in  (20)  can  be  used  to  evaluate  the  vv ,  but 
the  method  is  more  involved  than  evaluations  using  other  methods.  As  such,  it  is 
unlikely  that  this  method  can  be  incorporated  into  design  processes  involving  the 
selection  of  a  suitable  Q(jcj). 

2. 8. 2. 3  Fvaluation  of  a  Composite  Matrix  The  authors  of  (14)  offer  a 
necessary  and  sufficient  condition  for  the  stability  of  a  polynomial  matrix.  Changes 
in  the  presentation  as  made  by  the  authors  of  (14)  are  made  herein  to  make  the 
material  suitable  to  this  section. 

Define  P(,(ju>)  as  that  polynomial  matrix  resulting  from  the  product  of  Pe(jw) 
and  the  characteristic  Equation  (least  common  denominator  polynomial)  of  Pe(jw). 
Then,  for  any  nonsingular  l  x  /  polynomial  matrix  Pg(ju>),  a  diagonal  matrix  D ( jcd ) 
can  always  be  determined,  such  that  for  rmax  =  maximum  power  of  u>  in  the  nu¬ 
merator  terms  of  Pg(ju;), 


Pe(j")  =  B(iu){Ilxl(j^Ymax  +  ^m„-,(ju;)r— -1  +  ...  A,(ju,)  +  To]  (2.49) 
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where  the  maximum  power  of  the  polynomials  of  any  ith  column  of  P^jiu)  is  along 
the  ith  diagonal  of  D(jw>).  ft  follows  that 

de<(Pg(ju;))  =  det{D(ju)[Il>a(juj)TTnax  +  + ...  Ai(juj)  +  A0\)  (2.50) 

Now,  an  equivalent  expression  for  the  braketed  term  of  the  previous  equation 
is  given  by 

0/x/  •••  ®ixi  —  A0 

hxi  •••  0(x/  —  M 

0/X/  •••  hxl  ~Armai~.  I 

Hence,  it  is  concluded  that 

det( P'Jju;))  =  det(D(ju’))det((ju)Iix/  -  L)  (2.52) 

and  for  rmin  equal  to  the  smallest  power  of  w  in  D(ju/),  that  eigenvalues  of  Pg(jw’) 
are  given  by  the  eigenvalues  of  L  after  excluding  rmax  —  rmm  zero  eigenvalues  added 
by  the  formation  of  the  L  matrix.  The  determination  of  whether  the  vX]  are  MP  can 
be  made,  therefore,  by  determining  if  the  remaining  eigenvalues  of  L  lie  in  the  LHP. 

For  example,  for  the  following  Pe(ju/)  polynomial  matrix 


2-37 


the  det((juj)I  —  L )  yields  eigenvalues  of  0,  -1,  -2,  and  -  3.  Excluding  rmax  —  rTO,„  = 
2  —  1  zero  eigenvalues  (i.e.,  excluding  the  only  zero  eigenvalue),  we  determine  that 
eigenvalues  of  det(P'e(}uj))  are  -1,  -2,  and  -3  and  hence  the  vXJ  are  MP. 

Again,  the  method  appears  useful  in  screening  out  inappropriate  selections  of 
Q(juj),  hut  it  appears  unlikely  that  the  method  can  be  incorporated  into  the  selection 
process  of  Q(jcj). 

To  conclude,  methods  are  available  for  determining  if  the  (m,j  )  of  Pg^ju?)  are 
MP  by  evaluating  ^  ’?  ae<.;iJe(j_. )).  Il^wev-r,  si..v.e  the  available  methods  appear 
to  have  applicability  only  after  the  selection  of  the  Q(ju;)  compensator  is  made,  the 
methods  can  be  expected  to  see  little  use.  If  the  evaluation  is  done  after  the  selection 
of  the  Q(ju;)  matrix,  numerical  evaluations  of  det( Pe(ju>))  are  simpler  and  quicker. 
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2.9  Selection  o/Q(jw)  for  MIMO  QFT  Designs 


In  this  section,  QFT  design  requirements  are  examined  to  determine  optimal 
Q(jw).  As  the  intent  of  this  study  is  to  concentrate  on  the  selection  of  a  Q(ju;)  that 
minimizes  the  difference  between  the  two  sides  of  Equation  (1.1),  it  is  imperative 
that  P(jcj)  be  first  fixed  (  P(jiu)  €  Hence,  for  purposes  of  this  study  the 

objective  is  established  by  specifying  that  Q(ju>)  is  to  minimize  Equation  (1.1)  with 
respect  to  the  healthy  plant. 

The  first  criterion  for  Q(jw)  is  that  the  product  of  P(ju>)  and  Q(ju>)  be  a 
square  matrix.  This  constraint  is  stated  first  since  it  allows  us  to  proceed  under  the 
asumption  that  P(ju/)Q(ju>)  is  a  square  matrix.  Obviously,  for  any  l  x  m  matrix 
P(jw),  we  can  always  find  an  m  X  /  Q(ju>)  that  yields  a  square  /  X  l  Pe(ju;). 

Before  proceeding  to  the  second  criterion,  some  matrix  equivalents  are  estab¬ 
lished.  While  we  proceed  by  use  of  an  example,  similar  statements  can  be  made  for 
any  nonsquare  P(jcu)  and  Q(ju>)  matrices  whose  product  is  square  by  applying  the 
Binet-Cauchy  formula  (7:5-19). 

Consider  a  2  x  3  matrix  of  transfer  function  elements  P(ju>),  such  that 

/  \ 

.  /  Pn  P iz  P i3 

P(jw)  = 

\  P2l  P22  P23  ) 

and  a  generalized  compensator  Q(jw),  such  that 


/ 

\ 

?n 

<7l2 

Q(jw)  = 

921 

922 

K  931 

932  y 

The  product  P(jw)Q(ju;)  is  given  by 
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.  .  [  Pll<?ll  +  Pl2<721  +  Pl3<?31  Pi  1  <?12  +  Pl2?22  +  Pl3<?32  |  _  . 

P(jw)Q(jw)  =  (2.53) 

\  P21?ll  +  P22<721  +  P23?31  P21?12  +  P22<?22  +  P23?32  ) 


The  right  hand  side  of  Equation  (2.53)  is  rewritten  to  yield 


P«£?«  —  P1Q1  +  P2Q2  +  P3Q3  +  PaQa  +  P5Q5 

1=1 

-‘rP&Qe  +  E7Q7  +  P&Qs  +  P9Q9 


(2.54) 


Note  that  Equation  (2.54)  cannot  be  set  equal  to  the  right  hand  side  of  Equation 
(2.53)  without  first  defining  an  equivalence  relationship.  This  is  done  later.  For  now, 
continue  by  noting  that  the  P,  and  Q,  of  Equation  (2.54)  are  given  by 


PiQ  1  =  P(s) 


^  ?lll  ?ll2  ^ 

9l21  ?l22 

0  0 


Pi  1  P12 


\  / 


9lll  9ll2 


P21  P22  /  \  7l21  <7l22  / 


(2.55) 


P2Q2  =  P(s) 


72h  <72,2 

0  0 

V  <72,1  ?232  J 


( 


?2  n  <7212 

\  7231  72j2  J 


(2.56) 
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P3Q3  —  P{s ) 


0  0 
<?321  7322 

^  7331  <?332  ) 


'  Pl2  Pl3 
^  P22  P23 


732i  7322 

933 ,  7332  y 


'  <?4n  0 

(  0  q5l2  ' 

PaQa  =  P(s) 

0  <7422 

PsQs  =  P(s) 

<?521  0 

^  743,  <?432  ) 

^  95J]  7532 

/  \ 

(  \ 

76,,  <?6i2 

77,,  q 712 

P&Q6  =  P(s) 

0  <7622 

P7Q7  =  P(s) 

<7721  0 

^  ^631  0  y 

^  0  9732  J 

(2.57) 


( 

\ 

(  n 

\ 

78,, 

0 

0 

?9i2 

PsQs  =  P(s) 

782, 

?822 

P9Q9  =  P(s) 

7922 

0 

7832  1 

^  993] 

0 

/ 

Note  that  only  the  first  three  products  can  actually  be  formed  using  a  square 
Px.  The  2x2  P,  for  1  E  (4, 5, 6, 7, 8, 9)  do  not  truly  exist.  Let  the  3x2  Q,  matrices 
be  denoted  as  Q\.  Next  observe  that  each  of  the  Q[  for  i  E  (4, 5, 6,  7, 8, 9)  is  just  a 
linear  combination  of  the  Q't  for  i  E  (1,2,3).  As  such,  only  the  first  three  elements 
of  equation  (2.55)  are  needed  to  span  the  space  of  achievable  Q(ju>).  In  general 
for  an  /  x  m  P(jcj)  and  an  m  x  l  the  Binet-Cauchy  theorem  yields  m2  PXQX 

combinations  of  which  only  m  are  needed  to  span  the  space  of  all  possible  P(ju;)Q(ju;) 
products. 

That  only  the  first  three  PtQt  for  i  E  (1,2,3),  are  needed  to  span  the  space  of 
possible  P(ju;)Q(ju>)  is  verified  by  forming  the  sum  of  the  P,Q,  of  equations  (2.55) 
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through  (2.57).  The  sum  results  in  Equation  (2.54)  with  the  elements  of  Equation 
(2.54)  given  by 


9n 

-  [*?ln 

+ 

92, ,] 

<712 

=  [<?ll3 

,  + 

9212] 

<?21 

=  [?121 

.  + 

?321] 

<722 

=  [?12: 

,  + 

93m] 

<73 1 

II 

1  + 

<7331  ] 

<732 

=  (<723; 

i  + 

9332  ] 

As  such,  the  Pt  for  i  6  (4,  5, 6,  7, 8, 9)  are  eliminated  from  further  consideration. 

The  purpose  in  breaking  up  the  product  is  to  define  a  set  of  P,Q,  matrices  such 
that  any  Pe(jw)  =  P(jw)Q(ju;)  can  be  formed  as  a  linear  combination  of  several  P,Q, 
matrices.  The  objective  is  to  define  selection  criteria  for  Q(ju>)  on  the  on  the  basis 
of  the  square  PlQ,  matrices.  A  desirable  Q(ju>)  is  then  formed  using  the  desirable 
Q,.  Note  that  in  seeking  to  achieve  this  objective  that  Equation  (2.55)  potentially 
yields  m2  *  Pe(jcu),  since  Equation  (2.55)  contains  m2  versions  of  Pe(ju>).  If  only  the 
P,  for  i  €  (1,2,3)  are  retained,  then  the  relationship  between  equations  (2.54)  and 
(2.55)  is  given  by 

P(jw)Q(jw)  =  -  jr  P,Q„  n  =  3  (2.59) 

n  PP 


Note  however,  that  the  value  of  n  in  Equation  (2.59)  is  reduced  if  additional 
P,  are  eliminated.  For  example,  if  P2  were  not  invertible  or  failed  to  satisfy  the 
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remaining  QFT  criteria,  then  P2  would  be  eliminated  from  further  consideration 
and  the  value  of  n  in  Equation  (2.59)  reduced  to  n  =  2. 

The  second  criterion  is  that  the  product  P(ju.')Q(ju;)  =  Pe(ju>)  be  non-singular 
over  the  entire  range  of  plant  uncertainties.  A  sufficient  condition  for  implementing 
this  requirement  is  to  require  that  each  of  the  P,Q,  combinations  be  nonsingular. 

The  impact  of  a  P,  matrix  for  i  £  (1,2,3)  whose  determinant  det(P,)  does  go 
to  zero  for  certain  plant  parameter  variations  is  next  determined.  Note  that  this 
criterion  is  applied  to  functions  in  the  P,  matrices  and  not  to  evaluations  of  the 
functions  at  certain  frequencies. 

Assume  a  P,,  i  £  (1,2,3)  is  singular.  If  P,  is  singular,  then  for  some  plant  in 
the  set  of  possible  plants,  det(Pt)  =  0.  As  each  P,  is  square,  the  product  rule  for 
determinants  applies.  Hence,  for  P,Q,  combinations  where  i  £  (1,2,3), 

det{P,Q,)  =  det(Pi)dct(Ql)  (2.60) 

Based  on  the  above  result,  we  conclude  that  if  det(P,),  i  £  (1,2,3)  does  go  to 
zero,  then  no  Q,  eliminates  the  zero.  Conversely,  it  also  follows  that  if  these  three  P, 
are  nonsingular,  then  any  nonsingular  Q,  ensures  that  the  product  P,Q,  satisfies  the 
second  criterion.  To  proceed,  it  must  be  assumed  that  at  least  one  P,  of  Equation 
(2.59)  is  invertible. 

The  third  criterion  noted  in  Chapter  I  is  that  the  product  P(juz)Q(jiv)  satisfy 
a  diagonal  dominance  condition  as  w  -*  oo  over  the  entire  range  of  plant  parameter 
variations.  This  condition  ensures  that  the  diagonal  transfer  functions  dominate  the 
plant's  response  at  high  frequencies.  A  more  restrictive  version  of  this  requirement  is 
Rosenbrock's  dominance  condition  that  requires  that  the  diagonal  transfer  functions 
dominate  the  plant  s  response  at  all  frequencies. 


Note  that  diagonal  dominance  in  QFT  differs  from  that  usually  referred  to 
in  linear  algebra  courses.  In  this  case,  diagonal  dominance  is  that  characteristic  of 
a  matrix  that  is  observed  when  the  eigenvalues  of  the  matrix  satisfy  Gershgorin’s 
"circle  theorem"  (17:386-387). 

The  application  of  the  criterion  to  Pe(ja;)  or  to  the  Px  i  £  (1,2,3)  is  readily 
done  by  evaluating  transfer  function  elements  of  the  matrices  over  a  desired  frequency 
range  and  then  applying  magnitude  conditions  detailed  in  (7)  and  (3).  For  2  x  2 
matrices,  the  diagonal  dominance  condition  requires  that 


|PenPe22!  >  |Pe12Pe21|  (2.61) 

For  matrices  of  any  higher  dimension,  (7)  provides  the  necessary  information.  It 
follows  that  if  any  P,Q ,  that  is  diagonally  dominant  is  added  to  any  other  PjQj  that 
is  also  diagonally  dominant,  then  the  resulting  sum  is  also  diagonally  dominant. 

Once  a  Qx  matrix  is  selected,  then  the  matrix  product  PXQX  is  evaluated  over 
come  frequency  range  to  determine  the  extent  to  which  each  P,Qi  combination  is 
diagonally  dominant.  As  such,  we  are  able  to  select  for  QFT  design  purposes,  the 
diagonally  dominant  P,Q.. 

Ideally,  the  selection  process  for  Q(ju>)  yields  a  Pe(ju>)  that  satisfies  Rosen- 
brock’s  diagonal  dominance  condition.  In  such  cases,  design  efforts  are  reduced  and 
the  resultant  design  is  simplified.  However,  selecting  Q(jw)  to  primarily  accomplish 
this  objective  entails  abandoning  efforts  to  select  Q(jw)  on  the  basis  of  a  prespecified 
P,,(ju,’)  matrix. 

The  fourth  criterion  requires  that  roots  of  c/e<(Pe(jw))  lie  in  the  LHP  (i.e.,  that 
flie  =  rXJ  be  MP,  where  the  mXJ  are  elements  of  P'^jo,’)).  This  criterion 

cannot  be  applied  when  the  plant  transfer  matrix  is  nonsquare.  However,  in  the  form 
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of  a  set  of  squared-down  versions  PtQ ,,  the  criterion  is  readily  applied.  As  such,  the 
P,Qt  combinations  can  be  evaluated  to  weed  out  undesirable  choices  of  Q,. 

A  summary  point  about  the  four  QFT  criterion  is  now  made.  Of  the  four 
criteria,  three  are  GO-NO  GO  type  of  evaluations.  Only  one  of  the  four,  the  third 
criterion  concert  ing  diagonal  dominance,  may  be  said  to  be  qualitative  in  nature. 
As  such,  it  is  believed  that  the  selection  process  should  concentrate  on  the  third 
criterion,  reverting  back  to  the  three  remaining  criteria  only  to  the  extent  necessary 
to  ensure  compliance. 

Once  desirable  Q,  are  determined,  a  single  Q(ja.>)  is  then  formed  as  the  weighted 
sum  of  the  expanded  Q\.  Weighting  matrices  IT,  whose  sum  equals  7/x;  are  used  for 
this  purpose.  The  IT,  are  restricted  to  the  use  of  scalar  values  in  order  to  preclude 
further  changes  in  the  MP  characteristics  of  the  equivalent  plant.  The  QFT  design 
process  is  then  entered  using  a  single  Q(ju;)- 

The  fact  that  choosing  the  IT,  to  satisfy  IT,  =  /jx,  yields  P0(j.o)  is  readily 
demonstrated  for  the  case  of  invertible  P,  by  the  following  set  of  equalities.  The  O, 
satisfying  Equation  (1.1)  are  given  by  Qt  =  P~l  Pe  since  the  P,  must  be  invertible. 

P(Q\  IT,  +  Q'2W2  +  Q'W-3)  =  PQ\  VT,  +  PQ'2\V2  +  PQ' IT3 

=  P,  Pf1  P e IT,  +  P2P">  Pe\V2  +  P3P3-1  P'W  3 
=  Pe(VT,  +  W2  +  W3) 

=  pr 


Alternately  QET  designs  can  be  determined  for  each  of  the  P,Q,  for  retained  P, 
and  the  individual  designs  implemented  in  parallel.  With  the  QFT  process  completed 
u-ing  three  versions  of  f\Q,,  flesign  implementation  occurs  as  shown  in  Figure  2.3. 
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Figure  2.8.  The  QF'T  Compensated  Plant,  with  W  Weighting  Matrices 
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The  IT,  weightings  are  then  selected  to  achieve  an  overall  best  compensating  system. 
It  is  anticipated  that  if  tradeoffs  are  necessary  in  the  selection  of  the  VT,,  that  the 
tradeoffs  are  more  visible  in  a  post-QFT  environment. 

The  process  of  Hi  or  H^- norm  minimization  is  accomplished  in  the  absence  of 
other  design  criteria  for  establishing  tradeoffs  between  the  IT,.  As  each  P,  retained 
to  this  point  is  assumed  invertible,  it  follows  that  the  Qx  are  already  Hi  and  Hoc- 
norm  minimal  with  respect  to  the  healthy  plant.  As  such,  the  use  of  Hi  or  H^-norm 
minimization  is  used  to  address  the  entire  set  of  possible  plants.  As  there  is  no 
mathematical  formulation  available  for  this  purpose,  the  process  of  minimizing  the 
Hi  or  Hoc- norm  over  the  set  of  all  Pe(ju;)  G  Ve(\u>)  is  done  by  trial  and  error 
selections  of  the  IT,.  However,  the  process  is  applicable  only  if  transfer  function 
elements  of  the  retained  P,  are  stable  for  ail  plant  parameter  variations. 

Before  concluding  this  section,  the  idea  of  abandoning  a  prespecified  equivalent 
plant  and  instead  seeking  to  select  Q(ju?)  to  make  the  product  P(ju>)Q(ju?)  satisfy 
Rosenbrock’s  diagonal  dominance  condition  is  examined.  The  examination  is  con¬ 
ducted  by  examining  the  impact  of  two  different  Q(ju>)  on  the  diagonal  dominance 
characteristics  of  P(ju;)Q(ju;)  for  P,  G  (1,2,3). 

First,  a  diagonal  Q,  is  assumed.  For  explanation  simplicity,  it  is  assumed, 
without  loss  of  generality,  that  the  P,  being  used  is  P2.  The  product  of  Pi  and  an 
assumed  diagonal  Q2  is  given  by 


P7Q2 


Pn  P 12 
P21  P22 


72,,  0  ^ 
0  7222  ) 


(2.62) 


Carrying  out  the  product  and  applying  the  diagonal  dominance  condition 
liquation  (2.57),  the  following  is  obtained. 


2  17 


|(Pll<Z2„)(P22<7222)|  >  I(PX2?222)(P21?2u)I 
I<?2u?222|  i  P 1 1 P22  |  >  l?2,,?222|  |pi2P2l| 


and  it  is  concluded  that  a  diagonal  Q ,  provides  frequency  sensitive  scaling  of  each  side 
of  the  inequality,  but  does  nothing  to  alter  the  diagonal  dominance  characteristics 
of  the  P,  matrix. 

Now,  a  Qi  whose  numerator  terms  are  the  conjugates  of  numerators  of  elements 
of  Pi  and  whose  common  denominator  terms  are  the  same  common  denominator 
terms  of  P,  is  next  investigated.  To  express  both  the  P,  and  the  Q,  elements  in  a 
common  frame,  let  denp  represent  the  common  denominator  term  of  P,  and  denq  the 
common  denominator  terms  of  Q{.  Further  let  ntJ  represent  the  numerator  terms  of 
Pt-  Again,  for  explanation  simplicity,  it  is  assumed  that  the  Pt  being  used  is  P2. 

(  nn/denp  nl2/denp\  I  nTT/den,  rT^/denq\ 

P2Q2  =  _  _  (2.63) 

y  n2l/denp  n22/denp  j  y  rv^/denq  n^/denq  j 

Carrying  out  the  product  and  applying  the  diagonal  dominance  condition,  the 
allowing  is  obtained. 

W\\nun22n22  +  n\2n2in2i ^  1 2 1  ^  | 1 2 ri  1 1 P2 1^22  +  nuni2n22n2i  |  (2.64) 

where  the  common  denominator  terms  have  been  removed  since  they  only  alter  the 
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scaling  of  the  inequality.  But  as  the  right  hand  side  of  the  inequality  is  the  sum  of 
a  function  and  its  complex  congugate,  the  right  hand  side  satisfies 


1 2  *  /?e[rc12nnn21n22j|  =  |n12nnrc21n22  +  nun12n22n21|  (2.65) 


both  sides  of  the  inequality  are  equivalently  expressed  by 


I  | ^1 1  ^22 1 2  +  |rti2«2i|2  |  >  )  2  *  \j\nn\2  |nn|2  |n21|2  |n22|2  I  (2.66) 

Removing  the  excess  outer  absolute  value  bars  and  squaring  both  sides  of  the 
previous  inequality,  the  following  results 


^ i i ^22 14  +  2  *  | n i in22ni2n2i |2  +  |ni2n2i|4  >  2  *  |n12^nn2in22|2  (2.67) 


which  implies  that 


nnn22|4  +  |ui2n2i|‘<  >  0 


(2.68) 


and  it  is  concluded  that  this  Q,  guarantees  that  Rosenbrock’s  diagonal  dominance 
condition  is  at  least  satisfied  for  the  healthy  plant.  Hence,  we  are  able  to  define 
conditions  for  ineffective  and  effective  Q ,  for  improving  the  diagonal  dominance  of 
a  t\ .  However,  the  application  of  the  latter  Qi  entails  some  modification  since  a 
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selection  based  on  the  healthy  plant  may  drive  the  equivalent  plant  matrices  Pe(ju>)  £ 
7>e(jw)  to  be  NMP. 

2.10  Summary 

The  following  summarizes  material  in  this  chapter. 

1.  Elements  needed  to  determine  the  singular  value  decomposition  of  a  matrix  of 
transfer  functions  are  obtainable.  In  general,  however,  the  transfer  function 
elements  are  not  real-rational  functions  in  ju>  and  their  use  in  determining 
solutions  to  the  model-matching  problem  is  limited. 

2.  The  application  of  QFT  to  MIMO  plants  relies  on  defining  a  multiple-input, 
single-output  (MISO)  equivalence  mapping.  QFT  constraints  stem  from  this 
mapping. 

3.  Methods  are  available  for  determining  if  the  (m,y)  of  P"1^)  are  MP  for  a'l 
plant  parameter  variations  after  a  Q(ju;)  has  been  selected.  The  most  appro¬ 
priate  method,  however,  remains  the  numerical  evaluation  of  det( Pe(jo>))  for 
different  plants  Pe(ju;)  £  Ve(ju).  None  of  the  methods  available  for  determin¬ 
ing  if  the  (m,j)  of  P~A(ju;)  are  MP  is  readily  adaptable  to  the  forming  of  a 
mathematical  relationship  between  the  MP  criterion  of  QFT  and  a  desire  to 
minimize  either  the  H?  or  H^- norm  of  the  set  of  all  projected  plants  about 
some  nominal  plant. 

4.  Of  the  four  constraints,  only  diagonal  dominance  may  be  said  to  be  qualitative 
in  nature.  As  such,  it  is  believed  that  the  selection  process  should  concentrate 
on  satisfying  diagonal  dominance  over  the  largest  bandwidth  possible  (with 
the  objective  of  achieving  Rosenbrock’s  diagonal  dominance  condition),  and 
reverting  back  to  the  three  remaining  criteria  only  to  the  extent  necessary 
to  ensure  compliance.  To  do  so,  however,  requires  abandoning  a  prespecified 

Pe(jw). 
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5.  Breaking  up  the  product  of  P(ju)Q(ju;)  to  define  a  set  of  P,Qi  matrices  sat¬ 
isfying  Pe(ju>)  =  P(jw)Q(jw)  provides  a  number  of  insights.  QFT  designs 
can  be  accomplished  on  a  single  Q(ju>)  after  summing  the  nonsquare  Q\  or  on 
individual  PXQX  and  then  implementing  the  results  in  parallel. 

6.  As  only  invertible  P,  are  retained  for  QFT  design,  the  Q,  are  determined  using 
P,_1.  As  such,  each  Q ,  is  H2  and  //oo-norm  minimal  for  the  healthy  plant.  If 
the  design  objective  is  to  minimize  the  H2  or  //oc-norm  over  the  entire  set  of 
possible  plants,  then  it  is  suggested  that  the  process  of  minimizing  the  norm  be 
done  using  a  trial  and  error  approach.  This  process,  however,  effectively  rede¬ 
fines  the  plant  used  to  obtain  each  Q,  from  the  healthy  plant  to  that  nominal 
plant  which  minimizes  deviations  from  nominal  for  the  set  of  possible  plants 
with  respect  to  the  selected  norm  criteria.  Since  no  mathematical  formulation 
is  available  for  integrating  the  norm-minimization  requirement  with  a  sum¬ 
mation  of  P,Q ,•  equivalent  plants,  either  a  trial  and  error  aproach  in  selecting 
Q(ju>)  is  used  or  each  individual  Qi  of  PiQt  is  selected  to  norm-minimize  the 
P,QX  and  then  the  resultant  designs  implemented  in  parallel.  Finally,  to  apply 
norm-minimization  criteria,  it  is  noted  that  transfer  function  element  of  P(ju>) 
must  be  stable  for  all  plant  parameter  variations. 

7.  To  define  the  best  overall  compensating  system,  weighting  matrices  Wx  are 
introduced.  Elements  of  the  W,  are  chosen  to  satisfy  £  Wx  =  I  and  either 
overall  design  requirements  for  closed-loop  plant  performance  or  to  minimize 
the  norm-criteria  for  a  sum  of  P,Qi. 

8.  The  use  of  a  diagonal  Qi  with  an  invertible  P,  is  ineffective.  Diagonal  Q ,  do 
nothing  to  alter  the  diagonal  dominance  characteristics  of  a  P,  matrix. 

9.  The  use  of  a  Q,  whose  numerator  terms  are  the  conjugates  of  numerators  of 
elements  of  P,  guarantees  that  Rosenbrock’s  diagonal  dominance  condition  is 
satisfied  for  the  healthy  plant.  The  method,  however,  requires  modification  to 
guard  against  driving  equivalent  plant  matrices  Pe(ju>)  £  Ve{]ui)  to  be  NMP. 


III.  Flight  Controller  Design 


3.1  Examination  of  Plant  Characteristics 

In  this  chapter,  transfer  functions  are  expressed  using  Laplace  transforms.  Con- 
ditions  at  the  beginning  of  Chapter  II  are  satisfied;  the  frequency  domain  represen¬ 
tation  of  the  transfer  functions  is  given  by  substituting  jui  for  s. 

The  aircraft  plant  selected  for  compensation  is  that  of  the  AFTI/F16  in  level 
flight  at  0.9  Mach  and  at  an  altitude  of  20,000  feet.  The  objective  of  the  compen¬ 
sation  is  to  maintain  effective  plant  roll  rate  and  yaw  rate  response  despite  limited 
control  surface  failures.  The  linearized  plant  model  is  described  by  a  state  space 
equation  in  the  following  form 


x  =  Ax  4-  Bu 

(3.1) 

y  =  C x  -f  Du 

(3.2) 

where  the  positive  direction  of  roll,  yaw,  and  sideslip  is  as  shown  in  Figure  3-1.  A 
description  of  the  state  space  equation  variables  follows. 

6  is  the  aircraft  roll  angle  in  degrees 
,3  is  the  aircraft  sideslip  angle  in  degrees 
p  is  the  aircraft  roll  rate  in  degrees  per  second 
r  is  the  aircraft  yaw  rate  in  degrees  per  second 
is  the  aileron  deflection  in  degrees 
6r  is  the  rudder  deflection  in  degrees 
f>ht  is  the  horizontal  tail  deflection  in  degrees 
The  matrices  of  equations  3.1  and  3.2  follow. 


/? 

p 


6a 

i 

\ Skt 


u  = 


Sr 

\  6kt  ) 


0.0000 

0.0000 

1.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0345 

-0.3436 

0.0326 

-0.9976 

-0.0014 

0.0373 

0.0266 

0.0000 

-55.2526 

-2.8000 

0.1457 

-51.0500 

10.3950 

-50.7290 

0.0000 

7.2370  ■ 

-0.0232 

-0.3625 

-1.2501 

-5.8080 

-5.1371 

0.0000 

0.0000 

0.0000 

0.0000 

-20.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-20.000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

-20.0000 

1  0.0000 

0.0000 

0.0000  ^ 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

B  = 

0.0000 

0.0000 

0.0000 

20.0000 

0.0000 

0.0000 

0.0000 

20.0000 

0.00000 

^  0.0000 

0.0000 

20.0000  J 
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y  = 


c  = 


0.0000 

0.0000 


0.0000 

0.0000 


1.0000 

0.0000 


0.0000 

1.0000 


0.0000 

0.0000 


0.0000 

0.0000 


0.0000  N 
0.0000  J 


D  = 


/  0.0000 

v  0.0000 


0.0000  0.0000 

0.0000  0.0000 


The  matrix  of  transfer  functions  is  given  by  P  —  C(sl  —  A)~l  B  -f 
transfer  function  matrix  is  determined  to  be  of  the  form 


Pu  Pl2  Pl3 
y  P21  P22  P? 3  ) 


Letting  denp  be  the  characteristic  equation  of  P,  the  elements  of  P  are 

denp  =  (s  +  3.9083  ±;2.9614)(s  +  .0272)($  +  2.6973)(s  +20) 
Pu  =  (-1021)(s)(s  +  .3541  ±;2.9273)/den 
Pi 2  =  (207.91)(s)(s  -  4.6436)(s  +  5.0717)/den 
Pi 3  =  ( -1014. 58)(s)(s  +  .3749  ±  j'3.5777)/den 
P21  =  ( — 25)(s  +  .5461  1  _/3.2533)(s  +  1.1119)/den 
P22  =  (-1  16. 178)(.s  +  .3284  ±;.6928)(s  +  2.4822)/rfen 
P23  =  (—  102.742)(s  +  .6699  ±  >1.5474)(s  +  1.5372 )/den 


.  The 


(3.3) 
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For  the  nominal  (healthy)  plant,  all  eigenvalues  are  in  the  closed  LHP  and  the 
system  is  stable. 

When  plant  variations  occur,  for  example  as  a  result  of  control  surface  failures, 
the  plant  has  the  following  form.  In  this  alternate  form,  Aq  represents  the  effectiveness 
of  the  aileron  control  surfaces  6a,  k2  the  effectiveness  of  the  rudder  control  surface 
Sr,  and  k3  the  effectiveness  of  the  horizontal  tail  control  surfaces  Sht-  The  set  of  all 
possible  plants  is  represented  by  V  where  for  each  kt,  i  £  (1,2,3),  k,  £  [.8,  1]. 

p  _  (^i)Pu  {k2)pn  (k3)pi3 
\  (^i)P2i  {k2)p22  ( k3)p23 


The  squared  down  versions  of  V  are  next  examined  for  diagonal  dominance. 
Only  the  Px  for  i  £  (1,2,3)  need  be  examined.  This  is  done  to  establish  the  extent  to 
which  the  three  P,  are  diagonally  dominant.  The  examination  is  done  as  a  function  of 
frequency  as  plant  variations  occur.  The  condition  for  diagonal  dominance,  Equation 
2.62),  applies. 

Figures  3.1  through  3.3  result.  The  solid  line  plots  represent  the  extent  of 
diagonal  dominance  for  the  healthy  Pt.  A  positive  result  indicates  that  the  P,  is 
diagonally  dominant  at  the  indicated  frequency  while  a  negative  result  implies  the 
opposite  characteristic.  The  dotted  plot  of  each  figure  represents  the  minimal  amount 
of  dominance  present  as  the  k,  vary  between  1  and  .8  in  .1  step  increments. 

For  this  study,  few  points  are  needed  to  determine  diagonal  dominance  for  the 
given  plant  because  the  set  of  all  oossible  plants  described  by  Equation  (3.4)  satisfies 


3  ' 


dft 


(k,)pu  (kj)pij 

( k,  )p2,  ( kj  )p2j 


P'2 1  P  2j 


for  any  i  ±  j  and  i,  j  6  (1,2,3).  In  general,  a  much  more  complete  evaluation  of 
the  determinant  for  the  set  of  ail  possible  plants  is  required.  Note  that  it  is  the 
requirement  that  diagonal  dominance  be  satisfied  as  u  -*  x  that  drives  the  choice 
of  Pi  and  not  the  relative  area  above  the  line  denoting  where  the  dct(P, )  goes  to 
zero. 

Based  on  Figures  3.1  through  3.3,  it  is  concluded  that  Pi  and  P2  are  diagonally 
dominant  while  P3  is  not.  However,  the  transpose  of  P3  is  diagonally  dominant  and 
?!  is  theretore  retained  for  use  in  subsequent  design  steps. 

In  comparing  the  three  plots,  it  is  noted  that  the  points  at  which  all  three 
plots  go  to  zero  occur  at  the  same  frequencies.  Hence,  any  linear  combination  of  the 
three  plots  retains  the  same  frequency  points  at  which  |pq,p,22|  =  Ip.uPiai!  for  any 
Px.  This  is  an  important  point  relative  to  attempts  to  select  Q(ju?)  for  purposes  of 
ensuring  that  Rosenbrock’s  diagonal  dominance  condition  is  satisfied  over  the  set  of 
all  possible  plants.  The  current  situatior  indicates  that  selecting  a  Q,  that  satisfies 
Equation  (2.63)  results  is  a  P,Qt  that  still  retains  the  same  frequency  points  at  which 
IP...P.22I  -  |p.12P.'2i|- 

An  examination  of  the  determinants  of  each  of  these  three  P,  is  next  conducted. 
The  denominator  polynomials  c/ent,  den2,  and  den3  of  the  deteiminants  of  Vx ,  V2, 
and  V:i  respectively  are 


d<r,x  =  123815(s)(s  +  ,0272)(s  +  ,2901)(s  +  2.6973)(s  +  .3908  ±  j2.9614)(s  +  20 )4 
drn2  =  79.">35(s)(s  +  ,0272)(s  +  .2805)(s  +  2.6973)(s  +  .3908  ±  ;2.9614)(s  +  20)4 
dr  »3  =  1 39233(.s)(.s  +  ,0272)(s  +  .2932)(s  +  2.6973)(s  +  .3908  ±  j2.9614)(s  +  20)4 


As  each  of  the  denominator  polynomials  has  only  LHP  roots,  inverting  indi- 
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Figure  3.4.  Diagonal  Dominance  of  "PoCju;) 
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vidual  elements  of  the  P~l  yields  transfer  functions  that  are  MP.  For  the  particular 
plant  being  used  in  this  study,  the  denominator  polynomials  of  the  determinants  of 
Pi.  V2.  and  Vs  do  not  change  as  plant  variations  occur.  Hence,  if  the  healthy  plant 
P,  are  MP,  so  are  the  P,,  for  i  G  (1,2,3). 

3.2  Selection  of  Equivalent  Plant  Pe 

From  (4:23),  the  roll  mode  time  constant  is  to  satisfy  T,  <  l.Osec.  From 
(1:22),  a  requirement  that  ujnd  >  1  rad/sec  is  also  obtained.  These  requirements  are 
combined  with  the  maximum  allowable  oscillation  requirements  of  (4:78)  to  obtain 
desirable  functions  for  the  equivalent  plant  Pe. 

In  addition  to  the  performance  requirements,  the  difference  in  the  orders  of 
transfer  functions  in  Pe  are  selected  to  coincide  with  the  relative  order  of  transfer 
functions  in  P.  This  is  done  to  ensure  that  Q,  transfer  functions  are  strictly  proper. 

The  following  transfer  function  is  selected  for  tracking  command  inputs.  The 
function  has  a  damping  ratio  Cd  =  0.6  and  natural  frequency  u>„  =  7.5  rad/sec  for 
the  dominant  complex  pair  poles.  Hence,  peu  is  given  by 


peil  =  (562.5)/((s  +  10)(s2  +  9s  +  56.25)]  (3.5) 


The  following  transfer  function  is  selected  for  rejecting  undesirable  inputs.  The 
rejection  transfer  function  is  selected  to  satisfy  diagonal  dominance  requirements  of 
the  P,  while  retaining  the  denominator  polynomial  of  peu  . 


p,l3  =  (s  +  .001  )/[(s  +  I0)2(s2  +  9s  +  56.25)] 


(3.6) 
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Based  on  the  above  selections  for  peu  and  pei2,  the  desired  equivalent  plant  is 
formed  as 


Pe 


I  pen  Pe  12 

\  Pe  12  Pen 


(3.7) 


Table  3.1.  Figures  of  Merit  for  Pp 


Transfer  Function 

Peak  Value 

Final  Value 

tp 

t. 

Peil 

1.06142 

1.00000 

0.665 

0.894 

Pen 

0.118E-07 

0.267E-11 

0.211 

2.290 

Figures  of  merit  for  Pe  are  as  presented  in  Table  3.1.  The  settling  time  ts  and  the 
time  to  peak  tp  are  given  in  seconds.  Settling  times  given  are  for  when  the  time 
response  is  within  two  percent  of  its  final  value. 

The  time  responses  for  peu  and  pe,2  to  a  step  input  are  provided  as  Figures  3.4 
and  3.5. 

In  preparation  for  obtaining  a  QFT  design,  the  upper  and  lower  tracking 
boundary  functions  are  now  specified.  The  transfer  functions  chosen  for  the  up¬ 
per  tracking  function  Tu(s)  and  the  lower  transfer  function  establish  a  region 

of  acceptable  responses  that  satisfy  design  criteria.  In  addition,  to  satisfy  QFT  de¬ 
sign  steps,  the  relative  order  of  Tu(s)  must  differ  from  the  relath'**  order  of  Tl(s )  by 
at  least  2. 

I  he  following  function  is  established  as  the  upper  tracking  boundary  function 

P  {s). 
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Figure  Pel  2  Response  to  Command  Step  Input 


Figure  3.6.  Step  Input  Response  of  pei2 


Tu{s)  =  4(5  +  l6)/[s2  +  10.45  +  64]  =  4(5  +  16)/[s  +  5.2  ±  ;'6.0795] 


The  following  function  is  established  as  the  lower  tracking  boundary  function 

Tl(s). 

Tl(s)  =  6400/[s4  +  36s3  +  484s2  +  2880s  +  6400]  =  6400/[(s  +  8)2(s  +  10)2] 

The  frequency  and  step  input  time  responses  of  the  upper  and  lower  tracking 
functions  are  provided  as  Figures  3.6  and  3.7.  The  central  plot  in  each  of  these  two 
figures  is  that  of  the  desirable  equivalent  plant  pen . 

3. 3  Obtaining  and  Evaluating  the  2x2  Qj  Compensators 

The  Qi  compensators  are  next  obtained  using  the  healthy  2x2  versions  of 
each  Pi  for  i  €  (1,2,3).  For  the  given  plant,  each  of  the  P,  is  invertible  so  that 
Q ,  =  P~*  Pe.  Note  however, 

As  acceptable  Qt  were  not  able  to  be  determined,  this  section  continues  with  a 
description  of  steps  that  were  to  be  taken.  The  next  section  is  then  used  to  describe 
efforts  that  were  made  to  obtain  the  Qi. 

Assuming  the  Q,  are  obtained,  the  Q\  compensators  are  next  determined;  each 
Q[  is  of  the  form 

/  \ 

<7u  <?I2 

Qi  —  Q21  ?22  (3-^) 

^  931  932  } 

Having  determined  the  Q\  compensators,  the  product  pairs  PQ\  for  i  6  (1, 2,  3) 
are  formed  to  verify  achievement  of  Pe  for  each  cascaded  pair.  Alternately,  verifica- 
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tion  that  acceptable  Qi  are  achieved  is  available  by  obtaining  the  time  responses  of 
the  PtQ ,  pairs  to  step  inputs  and  comparing  the  step  responses  to  those  of  Pe. 

With  the  Q\  in  hand,  for  design  implementation  using  a  single  Q(ju)  matrix, 
the  selection  of  the  IF,  weighting  matrices  is  accomplished  and  Q(jcj)  is  formed  as 
the  weighted  sum  of  the  Q\.  The  IF,  are  selected  to  either  satisfy  design  requirements 
or  to  minimize  the  //2  or  //^-norm  over  the  set  of  all  possible  plants.  With  the  IT, 
weightings  determined,  the  QFT  design  is  then  accomplished. 

On  Obtaining  the  Q,  Compensators 

This  section  describes  efforts  made  to  obtain  the  Q ,  compensators  appropriate 
to  converting  the  healthy  plant  P\.  P2,  and  P3  into  Pe. 

In  (5)  the  following  operation  is  provided.  For  matrices  of  transfer  functions 
corresponding  to  the  state  space  realization  (.4,  B ,  C,  D),  define  [A,B,C,  D\  —  C  * 
[si  -  ,4]-1  *  B  +  D  =  P.  Then  P~l  is  given  by 

[.4,  B,  C,  D}~1  =  [A-  BD-xC,BD-\-D-xC,D~']  (3.9) 


Unfortunately,  the  above  equation  applies  only  where  D~l  exists.  As  D  for  a 
matrix  of  strictly  proper  transfer  functions  is  zero,  the  equation  is  not  appropriate 
for  obtaining  Q,  by  first  determining  P~ 1 . 

Because  of  the  unavailability  of  a  direct  solution  using  state  space  realizations, 
a  solution  method  based  on  determining  the  ratio  of  the  adjoint  of  P,  to  det(P, )  is 
attempted  ( 17:231-233).  The  method  is  applied  to  determine  each  Q,  as  the  product 
of  I\- 1  and  P,.  The  necessary  polynomial  product  manipulations  are  carried  out  on 
MATRIX*.  Results  of  the  manipulations  yield  what  appeared  to  be  acceptable 
in  verses. 
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The  Q ,  thus  determined  are  then  tested.  Results  indicate  that  the  Qx  do  not 
result  in  MP  [PxQi\~l-  The  inherent  difficulty  seems  to  lie  in  numerical  inaccuracies. 
While  some  of  the  inaccuracies  are  attributable  to  having  to  round  off  elements  of 
the  qXJ,  other  inaccuracies  are  determined  to  be  related  to  conversions  between  state 
space  and  transfer  functions  of  the  realizations  in  MATRIXx. 

For  example,  using  the  “eig”  command,  the  eigenvalues  of  the  A  matrix  given 
in  this  chapter  are  determined  to  be  the  values  indicated  in  column  1  of  the  following 
table.  Obtaining  first  the  transfer  function  form  of  the  state  space  realization,  and 
then  obtaining  the  roots  of  the  denominator  polynomial  (the  characteristic  function), 
the  roots  are  as  indicated  in  column  2.  Hence,  even  prior  to  manipulating  the  poly¬ 
nomials  to  form  the  adjoint  and  determinant  needed  to  obtain  the  P~: ,  a  noticeable 
error  is  present.  The  size  of  the  error  is  significant  in  that  its  size  is  the  same  as  the 
accuracy  of  the  linearized  plant  model.  It  is  believed  that  the  primary  reason  for  the 
failure  of  this  method  lies  in  such  numerical  precision  problems. 


Table  3.2.  Characteristic  Eigenvalues  of  P 


Eigenvalues  of  A 

Roots  of  characteristic 

-0.0272  -  jO.OOOO 

-0.0272  +  jO.OOOO 

-2.6973  +  jO.OOOO 

-2.6973  -  jO.OOOO 

-0.3908  -  j2.9614 

-0.3908  +  j2.96!4 

-0.3908  4-  J2.9614 

-0.3908  -  j2.96 14 

-20.0000  +  jO.OOOO 

-19.9998  -  jO.OOOO 

-20.0000  +  jO.OOOO 

-20  0001  +  jO.OOOl 

-20.0000  +  jO.OOOO 

-20.0001  -  jO.OOOl 

The  extent  of  the  impact  of  the  inaccuracies  is  provided  in  the  following  table. 
All  three  columns  of  the  table  are  roots  of  the  numerator  polynomial  of  <7in  whose 
largest  term  is  of  the  order  1  x  10H.  Column  entries  have  been  rounded  to  four 
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significant  digits.  Column  1  are  roots  determined  when  numerical  precision  is  kept 
at  16  significant  digits.  Column  2  are  roots  achieved  when  all  numerators  terms  with 
magnitude  less  than  1  x  10-6  are  rounded  to  zero  at  several  steps  in  the  computation 
process.  Column  3  are  roots  achieved  when  all  numerator  terms  with  magnitude 
less  than  1  x  10-3  are  rounded  to  zero.  The  wide  variability  in  numerator  roots 
determined  by  using  this  second  method  make  the  appropriateness  of  the  method 
questionable  in  this  case. 

The  third  method  attempted  is  to  solve  for  the  Q ,  directlv  by  forming  an 
equation  of  the  form  Ax  =  0  and  then  determining  the  solution  space  of  x.  This 
is  the  method  used  in  (12)  to  determine  the  Qt  matrices.  To  accomplish  this  an  A 
matrix  having  more  rows  than  columns  is  needed  so  that  x  is  known  to  lie  in  the 
nullspace  of  .4  (17:68).  To  this  end,  the  q,:  elements  of  Q,  are  defined  as  polynomials 
of  the  form 


ansn  +  an_]Sn  1  +  ...  ac 
i,J  bwsw  +  bw_lsw~l  +  ...  ba 


(3.10) 


under  the  assumption  that  the  a ,  and  6,  are  real  constants.  Both  sides  of  the  above 
equation  are  next  multiplied  by  the  denominator  polynomial.  The  resultant  left  hand 
term  is  then  substracted  from  both  sides  «,<>  yield 


[a„sn  +  an_i sn  -+-  a0]  ~  (jij[bwsw  -f  bw~isw  T  ...  b0]  —  0  (3.11) 


Evaluating  the  above  equality  at  a  minimum  of  n  +  w  +  2  frequencies,  q  > 
n  +  w  +  2,  the  resultant  set  of  linear  equations  is  written  in  matrix  form  as  follows. 
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Table  3.3.  Factors  of  q i,, 


No  Roundoff  Made 

10  6  Roundoff 

10  3  Roundoff 

0 

0 

0 

+6.154D-9,  -2.697 

+53.89,  -2.697 

+84.93,  -2.697 

-144.7,  -.2901 

-9.431,  -.2901 

-8.965,  -.2901 

-.02719,  -25.01 

-.02719,  -2.475 

-.02719,  -2.475 

-5.072,  +4.642 

-.3908  +  j2.961 

-.3908  +  J2.961 

-.3908  +  j2.961 

O 

00 

» 

to 

Id 

05 

-.3908  -  j‘2.961 

-.3908  -  j2.961 

-4.500  +  j6.000 

-4.409  +  j6.000 

-4.500  +  j6.000 

-4.500  -  j6.000 

-4.499  -  j6.000 

-4.500  -  j6.000 

-9.996  +  j. 07495 

-10.67  +  j 1 .37 1 

-9.753  +  J2.099 

-9.996  -  j. 07495 

-10.67  -  jl. 371 

-9.753  -  j2.099 

-14.31  +  j2.382 

-12.59  +  j5.381 

-11.08  +  J6.536 

-i 4.31  -  j2.382 
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To  obtain  the  q;j(juq)  requires  evaluating  the  product  P+(ju.’)Pe(ju;)  at  fre¬ 
quency  u.’,  and  extracting  the  ij  element  resulting  from  the  evaluation. 

To  obtain  a  set  of  vectors  spanning  the  nullspace  of  .4,  the  modified  Hermite 
normal  form  of  A  is  obtained.  The  “rref”  command  in  MATRIX*,  is  used  for  this 
purpose.  With  the  A  matrix  in  Hermite  noraml  form,  only  those  columns  whose 
diagonal  element  is  0  span  the  nullspace  of  .4.  A  “  -1  ”  is  inserted  at  the  diagonal 
position  to  yield  the  spanning  vector. 

In  uung  a  routine  based  on  this  third  method,  only  those  spanning  vectors 
whose  coefficients  are  real  can  be  used  (i.e.,  the  solutions  must  satisfy  the  assumption 
that  the  constants  of  equation  3.10  are  real).  Due  to  numerical  imprecision,  spanning 
vectors  have  typically  been  considered  real  if  the;,  imaginary  components  were  on 
the  order  of  10~10  or  less.  However,  in  running  the  routine,  imaginary  components 
were  always  obtained  with  magnitudes  well  above  the  typical  cut-off  (on  the  oidcr 
of  10-1  or  10~2).  As  such,  all  solutions  obtained  by  the  above  method  were  deemed 
unacceptable.  The  only  explanation  that  is  offered  to  explain  failure  of  this  algorithm 
is  that  the  product  P-1  (ju)P€(jcu)  was  ill-conditioned  in  the  frequency  region  of 
i nt crest. 

A  scluvion  using  a  curve-fitting  approximation  to  the  qij(juq)  points  determined 
above  was  not  attempted  due  to  time  considerations.  Curve  fitting  approximations 
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would  require  not  only  iteratively  fitting  the  curve,  but  also  establishing  t he  ac¬ 
ceptability  of  successive  approximations  using  the  QFT  criterion  on  the  product 

P(j-)Q(>’)- 

Finally,  a  solution  based  on  determining  the  P~l  by  first  augmenting  the  nu¬ 
merator  terms  of  transfer  function  elements  of  P,  with  sufficient  zeros  at  infinity 
(e.g..  by  multiplying  each  p%1  by  (s  4-  106)2)  to  make  the  D  matrices  of  each  P,  in¬ 
vertible  so  that  Equation  (3.9)  ean  be  applied  was  also  not  attempted  due  to  time 
considerations. 
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IV.  Conclusions  and  Recommendations 


This  thesis  focused  on  the  selection  of  a  frequency  sensitive  weighting  matrix 
Q(k)  to  convert  an  original  non-square  /  x  m  plant  transfer  matrix  P(jiu)  into  a 
square  plant  matrix  appropriate  for  multiple-input,  multiple-output  (MIMO)  quan¬ 
titative  feedback  technique  (QFT)  design  work.  The  plant  was  assumed  to  have 
more  inputs  than  observable  outputs.  The  method  of  pre-specifying  an  /  x  /  desir¬ 
able  plant  transfer  matrix  Pfl(ju.')  for  determining  the  necessary  nonsquare,  rn  x  / 
weighting  matrix  Q(j^)  was  investigated. 

W  hile  a  number  of  conclusions  concerning  the  selection  of  an  optimal  Q(j^') 
were  established  and  summarized  at  the  end  of  Chapter  II.  the  basic  problem  of 
integrating  QFT  design  criteria  with  norm-minimization  criterion  remains. 

Difficulties  encountered  in  obtaining  Q,  indicate  that  it  may  be  more  advisable 
to  pursue  optimization  methods  that  yield  a  Q(ju>)  directly.  The  processes  followed 
in  this  study  to  obtain  Q,  =  P~x  Pe  did  not  provide  suitable  Q,  for  the  given  plant 
and  Pe(jw’);  numerical  precision  difficulties  encountered  in  this  study  were  sufficient 
to  introduce  RHP  factors  into  c/e£(Pe(ju;))  over  the  set  of  plant  parameter  varia¬ 
tions  even  though  det( P(ja/))  had  no  RHP  factors  over  the  set  of  all  possible  plant 
parameter  variations. 

The  following  recommendations  for  areas  of  future  study  are  made: 

1.  Investigate  the  use  of  an  Hi  or  H ^  optimal  controller  Q(jw)  that  robustly 
stabilizes  all  plants  P(jcu)  in  a  set  of  possible  plants  P(j u)  for  the  case  where 
the  disturbance  d2(jcv)  is  zero. 

2.  Investigate  the  use  of  a  Q,  whose  numerator  terms  are  the  conjugates  of  some 
element  Px  £  V,,  the  set  of  squared-down  versions  of  that  plant.  Determine 
which  Pi  G  Vi,  if  any,  allow  a  cascaded  PtQ,  pair  to  either  reduce  the  band- 


4-1 


width  over  which  diagonal  dominance  is  not  achieved  or,  if  possible,  satisfy 
Rosenbrock's  diagonal  dominance  condition. 

3.  Continue  the  investigation  effort  of  this  study  after  identifying  an  acceptable 
method  to  obtain  either  the  P^1  or  the  Qx  matrices.  Determine  if  the  use  of 
the  pre-specified  Pe(ju,>)  improves  the  effectiveness  of  a  QFT  design  involving 
the  given  plant. 
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